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POREWORD 



This 1 b a combined report for the two projects at the Massachusetts Institute of PART I 

Technology which are sponsored by the Office of Naval Research under Contract N5orj60. 


Project on Machine Methods of Computation and Numerical Analysis 

Th'.a Project la an outgrowth of the activities of the Institute Committee on Machine 
Methods of Computation, established In November 1950. The purpose or the Project Is (l) to 
Integrate the efforts of all the departments and groups at M.I.T. w,‘o are wording with modern 
computing machines and their applications, and (2) to train men In the use of these machines 
for computation and numerical analysis. 

People from several departments of the Institute are taking part in the project. 

In the Appendix will be found a list of the personnel active In this program. 


Project Whirlwind 

This Project makes use of the facilities of the Digital Computer Laboratory. The 
principal objective of the Project la the application of an electronic digital computer of 
arge capacity and very high speed (Whirlwind l) to problems In mathematics, science, engi¬ 
neering, simulation, and control. 


The Whirlwind I Computer 

Whirlwind I Is of the high-speed electronic digital type. In which quantities are 
represented as discrete numbers, and complex problems are solved by the repeated use of 
fundamental arithmetic and logical (l.e., control or selection) operations. Computations 
are executed by Tractlona1-mlcrosecond pulses In electronic circuits, of which the princi¬ 
pal ones are (l) the flip-flop, a circuit containing two vacuum tubes so connected that one 
tube or the other Is conducting, but not both; (2) the gate or coincidence circuit; (3) the 
magnetic-core memory. In which blnax-y digits are stored as one or two directions of magnetic 
flux within ferro-magnetlc cores. 

Whirlwind I uses numbers of 16 binary digits (equivalent to about 5 decimal digits). 
This length was selected to limit the machine to a oraotlcal size, but It permits the computa¬ 
tion of many simulation problems. Calculations requiring greater number length are handled by 
the use of multiple-length numbers. Rapid-access magnetic-core memory has a capacity of 
32.768 binary digits. Present speed of the computer Is AO.000 single-address operations per 
second, equivalent to about 20,000 multiplications per second. 


Machine Methods of Computation and Numerical Analysis 


1. GbNKHAL COMMENTS 


Since the last report, arrangements have been made with the International 
Business Machines Corp. (IBM) for the Installation of a type 704 machine at the Institute 
In the spring of 19S7 and for the organization of the educational and general scientific 
use of computing equipment at MIT Into a Computation Center, with quarters In the new 
Compton Laboratory, Heretofore, Vhlrlwlnd I has provided most of the computational sup¬ 
port for this program. Whirlwind will eventually be turned back to Lincoln Laboratory 
about July 1, 1957, and the 704 and other equipment In the Computation Center will make 
possible an expanded program of research and training In the use of electronic computers. 

This expanded program will Involve the participation of the other colleges and 
universities in New England, and a Planning Committee, with representatives from 23 In¬ 
stitutions, has been set up to work out details of the Joint program. IBM has provided 
funds for the appointment of part-time research assistants and associates, from MIT and 
also frem the other participating colleges, to assist In carrying out the research. It 
is hoped that the National Science Foundation (NSF) will provide funds for full-time re¬ 
search staff. 

The coming academic year, therefore, will be a transitional one, changing from 
the ONR-supported project which has been reported in the previous Quarterly Reports, to 
the new program, supported by IBM and NSF and participated In by other New Ingland ool- 
leges. The first change Is the cessation of the part called Project on Machine Methods 
of Computation; by September 1956 all 0NR part-time research aaslstantshlps will be ter¬ 
minated and the new IBM awards will take their place. The Digital Computer Laboratory 
and the use of Whirlwind I will continue during the next academic year, but the staff of 
Project Whirlwind will be gradually transferred to the payroll of the Computation Center 
so that, by July 1, 1957, ONR-supported work will have ceased. 

The scientific research, reported In this and the preceding Reports, will not 
cease, however. It will continue at an Increased pace, and the only difference will be 
In the source of financial support. In line with this continuity of research Into the 
utilization of computing equipment In science and engineering, we propose to continue 
reporting the activities of both the older Project Whirlwind and the newer Computation 
Center In the one quarterly Report. This Issue will have the same form as before, with 
reports from the Project on Machine Methods of Computation and from Project Whirlwind. 
Next quarter, and the succeeding four quarters, will have reports from Project Whirlwind 
and also from the MIT Computation Center; after that the reports will be solely from the 
Computation Center. Distribution of the Reports will, of course, change as the now 
activities develop. 

Philip M. Morso 

Director, MIT Computation Center 


A 
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2.2 Progress Reports 

DETERMINISTIC PROCESSES AND STOCHASTIC PROCESSES 

In this report we summarize work that has been in progress during most of the 
year concerning deterministic processes and probability. This is an area of research 
that begins in statistical mechanics and is dirseted toward the interpretation and use of 
probability theory. The results could be useful for computation techniques because of 
what they imply about pseudo-random numbers and Monte Carlo techniques. The details ars 
contained in definition and theorem form and will be published elsewhere. 

A dynamics in the physical sense can be represented mathematically by a one- 
narameter group of transformations 

{V, < <*.} 

defined on and to a space of points A . if a point w in XI represents the state of a 
physical system at time t=0, then T w represents the state of the system at time t. In 
statistical mechanics we take Cl to be phase space and T to be the Hamiltonian dynamics. 
Liouville's theorem says that any set of points in P are transformed by T in such a way 
that the volume measure of the set in phase space is preserved in time. Such s property 
for T is given the mathematical name "measure-preserving. ' For definiteness we will call 
a one-parameter group of measure-preserving transformations a deterministic process. The 
knowledge of the state of a physical system at any fixed time together with the 
knowledge of a deterministic process operating on the system will determine all future 
and past states. 

Contrast the idea of a deterministic process with the idea of a sequence of 
independent and identically distributed random variables. For simplicity let a random 
variable correspond to each integer point in time. There is no apparent dynamics that 
links the various sample values of these random variables. On the other hand, we will 
show that a deterministic process with special properties can be Incorporated into the 
definition of random variables so that almost no generality is lost and yet a determi¬ 
nistic interpretation can be given. Both probability and statistical independence be¬ 
come asymptotic properties of deterministic processes. We thus conclude that abstract 
probability theory may be useful solely because of the relative orders of magnitude of 
things measured and things measurable and because of the order of magnitude of the time 
between measurements relative to characteristic parameters of a deterministic process. 

In order to ^ain more insight into these ideas it is necessary to use two 

basic theorems in statistical mechanics that were postulated by physicists and later 
proved (in modified form) by mathematicians. Tbe mathematical condition that it has 
been found necessary to Impose on a group of transformations T 1 in order that both of 
these (modified) theorems hold is called "mixing". It is a group of transformations 
satisfying this condition that we wish to Incorporate into a definition of random varia¬ 
bles. The two theorems, themselves, are designated ergodic" and mixing" and are 
roughly described as follows. In the ergodic case, we consider a single trajectory or 
"orbit” of a point in phase space obtained by operating on this point by a group of 
measure-preserving transformations. The ergodic theorem asserts that if a measurable set 
of points in phase space is specified then all trajectories except those beginning in a 
set of measure zero ultimately enter the measurable set infinitely many times and spend, 
all told, an average amount of time there that is equal to the measure of the set. In 
the mixing case, we must consider a collection of trajectories. The weak mixing theorem 
states that for any measurable collection of trajectories the sete of points that the 
trajectories occupy at time t=0 and at large time t are such that the measure of t;heir 
intersection is approximately equal to the product of their measures. This result re¬ 
quires that large t is chosen outside of a temporal set of density zero. We see in the 
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first of these two theorems the gera of an Interpretation for probability, for proba¬ 
bility can be defined as the measure of a set in phase space. In the second of these 
theorems we see an interpretation for statistical Independence. 

Any definition of randosi variables requires a basic space of events XI , a 
Borol field, B, of sets on this space and a probability measure P defined for all sets 
in B. We may choose for these three entities the pnase space, the measurable sets in 
phase space, and the Lebesque measure. We can show that there exists u countably dense 
field F of sets in B in the sense that the number of sets in P is countable and every 
set in B can be approximated arbitrarily closely (symmetric differences being small) by 
sets in F. Sinco P is countable, the ergodlc theorem, stated with the some exceptional 
set of measure zero, applies to all sets in F. Thus if the dynamics T satisfies the 
m»xin r ; condition so that the ergodlc theorem is true, there oxists a trajectory that 
wnters all sets of F and spends in each set an average amount of time that is equal to 
the measure of the set. We can then consider the intersection of this trajectory with 
phase space. For each set in F there is a set within the trajectory and we can arbi¬ 
trarily transfer the measure from one to the other. With this device we can then view 
the trajectory (p as a new space of points, the sets F' so manufactured as composing a 
field of sets within <p , and the new measure P' as defined on F' . If necessary we can 
generate the smallest Hor vi field bF' containing P’and extend the measure P' to this 
field, so that the three entitles are available for defining random variables. 

Instead of defining a random variable on the measure space ( X} , B, P) as is 
customary, we define it on the measurable group ( , BF' , P'). We speak of a measurable 

group because (£ is now isomorphic to the transformation group representing the dynamics. 

A random variable X is simply defined as a measurable function on <p , distinct values of 
X mapping Inversely into sets of BF’. If we write 

<t> = {t* , -"t^l 

then for any point T 1 in th» sample space <P the random variable X takes on the value 
XCT 1 ). If we define the T-transiate of X by requiring that X^CT*) = X(T t+T ) then it can 
be seen that the mixing property implies for a wide class of randon variables that X- 
and X become statistically Independent with Increasing ?. 

Ws conoluds with s brief theorem suggestive of the type of theorem that can now 
be obtained. 

THBORKN 

Let X and Y be random variables dsflned on the measurable group (<£, BP', P') 
such that the smallest Borol fields Induced by X and Y in $ belong to F’ . Assume that 
F' Is Invariant under all rational transformations 7*6 <j> , r rational. 


2) For any S ? 0 and any real Borel sets B t , B 2 there exists a rational r such that 
|P’[X < Bj, Y r e B a ] - P'[X d Bj P'[Y r * Bj | C f 

Bayard Rankin 

A CALCULATION OF THE ENERGY BANDS OF THE GRAPHITE CRYSTAL BY MEANS OF THE TIGHT- 
BINDING METHOD 

The computation of energy band! in a two-dimensional model of graphite by 
means of the tight-binding method has been completed. Current effort is being devoted 
to the thesis write-up of this work. Thus the following report la only a brief resume' 
of the basis of the calculation and the major results which were obtained. 

The graphite energy band calculations were carried out assuming a structure of 
an infinite two-dimensional hexagonal lattice of Carbon atoms. The atomic orbitals used 
to make up the one-electron Bloch waves were the Hartree-Fock atomic orbitals of the 
ground state neutral Carbon 3 P configuration calculated by Jucya 1 . The Is and 2p orbi¬ 
tals were each fitted by a linear combination of three corresponding analytic Slater a- 
tomic orbitals, and the noded 2s orbital was fitted by a sum of three Is and three 2s 
Slater atomic orbitals 2 . The fits were of high accuracy and were made for ease in the 
calculation of the many integrals. In the sense that the energy band calculation made 
was a one-electron approximation, it was necessary to choose an effective potential for 
the crystal. The crystal potential was taken to be a super-position of spherically 
symmetric atomic potentials. The atomic potential in turn was assumed to be the Coulom- 
bic potential arising from the effective nuclear charge function, Z (r), centered on the 
atomic nucleus. The Z p function for the 3 P configuration has been P calculated by 
Freeman 3 , and it was because of this convenience that the particular configuration of 
Carbon chosen was used since it was not felt that the current knowledge of energy band 
solutions gave any other conclusive choice. For computational convenience, the Zp 
function was also fitted, a linear combination of four exponentials bel.ig used. 

All the one-electron, two-center integrals (overlap, kinetic energy and poten¬ 
tial) were done by the use of the usual |.rolate spheroidal coordinate analytic Integra¬ 
tion techniques suitable for Slater atomic orbitals 4 . Thus these integrals required the 
evaluation of auxiliary functions 3 * 6 followed by combinations of many terms. All the 
one-electron, three-center potential integrals (and also the two-center potential inte¬ 
grals) were evaluated by the spherical coordinate expanslon-about-another-center tech 
nlque, similar to that used by Lowdin 7 , and Barnett and Coulson®, which was described 
extensively in a previous report 9 . 


1) For any real Borol set Bj.any rational r, and for all t 

' h * re otherwise , 


It is perhaps illuminating to give some of the computational perspective in¬ 
volved in the over-all calculation. Inasmuch as the magnitude of numerical work in¬ 
volved was such that accuracy would have been difficult to maintain in any hand calcula¬ 
tion, most of the computational work was done on the high-speed electronic computer 
Whirlwind I. The nearly total mechanization of the problem, although eliminating almost 
entirely any possible random mistakes, had the disadvantage of tending to obscure possi¬ 
ble systematic mistakes. The latter shift of emphasis is one of the principle reasons 
that make the programming of a computer a non-trlvlal affair. A consequence is that 
logical simplicity of the computational procedure becomes a goal which is often in 
opposition to computational efficiency. There ie also a great deal of difficulty in de¬ 
vising adequate test procedures for computer programs, since cf neceeslty they must be 
tailored to the program itself. The ability of the programmer to cope with these computer 
problems develops mostly with experience. A large fraction of the time spent on the 
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present calculation was thus used learning how to obtain the full potentiality of a high¬ 
speed computer. 


In carrying out the calculation, the work fell Into stages for each of which 
special computer programs were written. These were: a program for the semi-automatic 
fitting of the atomic orbitals and the Z function^; a program for generating the two- 
center integral auxiliary functions and then automatically combining terms to give the 
Integrals between Kartree-Pock orbitals 5 ; a program for generating the atomic orbital ex¬ 
pansion functions necessary for the three-center potential Integrals 9 ; a program for per¬ 
forming the basic numerical quadratures of the three-center potential integrals 9 ; and a 
requantization and summation program for forming the appropriate three-center Integrals 
from the basic numerical quadratures 9 . The foregoing computer programs were sufficient 
machinery to prepare the basic two-center Hamiltonian and overlap Integrals which served 
as Input for the final master program which performed the energy band calculation. Ex¬ 
plicitly, the basic two-center integrals were of the form 

“ (f) ( 1 ) 


where the Ijf are atomic orbitals. 


is the neighbor-vector, and 


M(r) 



1 


- 1/2 V 2 



<*>- 

n^: 


(Overlap) 


(Hamiltonian) 


( 2 ) 


the sum In the potential term being over-all neighbor-vectors. In the preaent calcula¬ 
tion, the three-center potential terma were neglected when |£ b | or I ? b -f a I exceeded 
the fourth neighbor-distance In the sum over P b . 


The operation of the master energy band program then proceeded as follows. For 
a given value of the reduced wave vector, k, the program coaiputed the Hamiltonian and 
overlap matrix elements arising from the Bloch waves constructed from the atomic orbitals 
These matrix elements were made real by taking judicious linear combinations of the Bloch' 

waves as basis states, and had the form 


V - f®°« v 

\j (t) “ * VV < 

f J l»in (k • P > 

a \ * 


(3) 


As has been implied, the secular equation which was considered had as basis 
states linear combinations of the ten one-electron Bloch waves formed from the Is, 2s, 
and 2p atomic orbitals on each of the two atomic sites In the spatial unit cell. Be¬ 
cause of the reflection symmetry In the Hamiltonian operator of the two-dimensional 
graphite lattice, the secular equation one obtains from these ten states Immediately 
factors Into two independent equations, one of order eight arising from reflection sym¬ 
metric Bloch waves (sigma states), and the other of order two from reflection antisym¬ 
metric Bloch waves (pi states). Thus the master computer program was arranged to Inde¬ 
pendently calculate the energy band solutions arising from each of these secular equa¬ 
tions, but the final results were graphically superimposed. 

The major physical significance in a two-dimensional energy band calculation 
of graphite is the size of the gap between the five lowest (occupied) and the higher 
(excited) sigma bands, and in particular whether or not the two lower pi bands (valance 
and conduction), which are degenerate at one value of the wave vector, have their point 
of degeneracy in the sigma band gap. If such a sigma band gap Is large enough to In¬ 
clude all of the pl bands, then It follows that a reasonable approximation for computing 
the conduction properties of a three-dimensional graphite would be to Ignore the sigma 
states; (in the three-dimensional crystal the sigma and pl are no longer good” symmetry 
designations since the sigma and pi Bloch waves of alternate graphite layers interact, 
but the terminology Is still used.) Starting from the original work of Wallace 10 , this 
approximation has been the basis of all graphite energy band calculations with the ex¬ 
ception of that of Lomer who has also considered a two-dimensional graphite with sigma 
bands 1 . As will be seen later, the numerical work of Lomer Is believed to have a very 
serious approximation. 

The energy band solution of the present calculation is shown in fig. 2 where 
the horizontal dimensions are the edges of the basic triangle of the first Brlllouin 
zone as shown in fig. 1. Energy values of Interior points of the triangle were also 
found, but for brevity these values are not given here, since it was found that the in¬ 
terior energy bands vary smoothly from the edge band values. Furthermore the two lowest 
sigma bands arising almost entirely from the Is Bloch waves, are omitted from fig. 2 
since the bands are nearly independent of the wave vector and at an average value of 
-15.75 Rydbergs. It io observed that the highest occupied point In the sigma bands oc¬ 
curs at pt. 0, (lc=0), with a value of +.430 Rydbergs. The lowest excited sigma band has 
roughly a constant minimum value for wave vector values forming an approximate circle 
about pt. 0, the minimum value being about +1.158 Rydbergs. The degeneracy point of the 
pi bands falls at +.570 Rydbergs, a value representing the Formi level for zero tempera¬ 
ture. Thus It is clear that these results support to some extent the usual approxima¬ 
tion of neglecting the sigma states in calculations of graphite conduction properties. 
Moreover this calculation shows by virture of the overlapping of the sigma and pi bands 
that for a cohesive energy calculation, which depends on energy values for all wave 
vectors of the Brlllouis zone, the sigma states must be included. 


wher. th. terns In which t exceeded the ninth neighbor-distance were neglected Die 
progran naxt aolred the usual variatlonally-derlved secular aquation of tha form 

X V {)V jk (f) “ *1<*I 8 lj( C > V jk (C > (4 

J J 


and atorad for l.tar u.e tha elg.nT.luaa K (x). A naw value of tha .... rector «a. th 
lected and the g.naratlon and solution ol the secular aquation r.peatad until tha nr 
t values of thw wave vector were exhauataH iMn.ii. ___ v 


sel 


swuwi«i.auii m.11 U Buiui lull or l( 

set values of tba ears rector were exhausted. Plnally, for conrenlince; the progran 
displayed graphically on a photographic oscilloscope cross-sectlonal view. of the ener 
bands, *l(*> rs. t, for raluas of the .are rector along the edges of a basic non¬ 
repeating 30° -60° -90° trlangl# of tha fir.t Brlllouin none 


Consideration can also be made of the potential function used in the present 
calculation. The potential, which was formed by superposition of the Z function, 
clearly omits exchange effects. Slater has given a procedure for Introducing an approxi¬ 
mate exchange potential correction within the framework of the one-electron approxima¬ 
tion 12 . This exchange potential correction, If applied, would bo proportional to the 
cube root of the charge density of the occupied crystal wave functions. Consequently, 
one would expect the occupied bands to be lowered more In energy than the unoccupied 
bands since the unoccupied wave functions are orthogonal to the occupied wave func¬ 
tions . Thus it is plausible that a more careful consideration of exchange effects In 
the present calculation would only broaden the signs band energy gap and would leave the 
present results qualitatively unchanged. 
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One of the more striking features of the present results Is the smoothness of 
the energy bands. In fact those bands are similar to those obtained from the Slater and 
Roster interpolation procedure wherein the tight-binding method is used with neglect of 
all but a few nearest-neighbor integrals^. It is therefore of interest to oxamine the 
stability of the present results in view of the several possible simplifying approxima¬ 
tion* . 


Omission of the Is Bloch waves was found to have an over-all lowering and 
warping effect on the sigma bands in such a way that the sigma band gap was roughly re¬ 
duced by half. Thus the often-recognised importance of orthogonality is again emphasised. 

A second possible approximation made was the omission of all the three-center 
potential integrals. This left the pi bands nearly the same but made a very pronounced 
change in the sigma bands, again closing the sigma band gap down to about half but also 
lowering the gap so far that the pi band Fermi level no longer was included. 

Finally the solutions were examined with respect to the approximation of 
omitting the higher neighbor integrals. The stable solution, which Included up to ninth 
neighbor Integrals, was found to be only slightly warped when only up to fifth neighbor 
integrals were included in the sigma bands and when only up to third neighbor integrals 
were Included in the pi bands. However, further truncation of the neighbor integrals of 
the sigma bands caused violent changes. (It is for thlB reason that the calculation made 
by Lamer which Included only first neighbor integrals, is not believed to be valid.) In 
addition, the effects of truncating only Hamiltonian or only overlap neighbor integrals 
were investigated. It was found that the solution was sensitive to both Hamiltonian and 
overlap neighbor truncation to roughly the some extent and that the two effects were 
essentially additive. 

Thus the results of all the approximation tests Indicate clearly that the 
tight-binding method when used in a non-emplrlcal way must be carried out with con¬ 
siderable mathematical rigor in order that a meaningful solution will be obtained. 
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THE GROWTH OF FATIGUE CRACK3 

An oxact solution to tho problem of determining the rate of propagation of a 
crack under cyclic loading is exceedingly difficult, and analytical solutions may be ob¬ 
tained only by approximating the exact formulation, even for simple loading modes and 
specimen configurations. The specific caso chosen for study here Is the torsion of a 
long cylindrical bar of rectangular cross section when the crack is advancing on a plane 
through the axis of the cylinder. The following assumptions are made: 

1) The material is perfectly plastic, homogeneous, and isotropic; 

2) The bar has been strained to an extent that plastic flow has occurred 
over virtually tho entire cross section; 

3) Tho crack propagates in tho following Banner: 

a) The bar is twisted until the material a distance, Xs, in front 
of the crack reaches the fracture strain, 

b) The crack length, f* , then increases by the distance df . 

c) This advance relieves stresses and the angle of twist, 9, must 
be increased by d0 before the fracture condition is again obtained. 

Within the limitations of theso assumptions, the rate of crack growth, 0, when 
the crack is of length X, may be expressed as the Volterra integral equation, 

• 2 „ , X. * 
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X fl 1 8 ^ [ _1_ 

- 1 + (^T7 + I7 1 J 1 + 

8 8 8 n 1 


when f(l) goes from 1 to 0. Terms of the order X 2 are small compared with one and have 
been discarded. A trial solution of the form 

inserted with the equation and iterated once to evaluate the constant C yields the re¬ 
sult that C has a value approximately between two and four. Thus we can expect j(X) to 
decay to 96% of stepdy state value after the crack has progressed about 8X a to 16X S . 

To evaluate the integral equation numerically, the oquatlon is written in 
finite difference form with spacing, h, where 


If the integral is approximated by a straight line between intervals, the equation takes 
the form 

#<nh) = l-0 8g * “ss ( ,lW ' X j] + X s A(n,v)0(vh) 


A(n,0) = 2 K(n ,0), 

(A(n,v) = h K(n,v) 


where the variables have been non-dimensionallzed. 

A steady-state solution to this equation may be obtained. 


0 «* _ v 2 „ 2 X 
ss In —(1 —) 


Using this result, tho problem may be re-fomuilated as 


0(X) - 0 




x . / [ 1 + TT- (X+ V] 2 ] 


A(n,n) = - K(n,n) , 


K(1 ,■)•! + 


/_i_y 

Uh-abfX I 


When the summation is evaluated at v=0 and v=», the difference equation becomes for n>l, 


i<“> ■ Sr|(s^ - *-"(=5^;)) [' - 


(“">.) J * 2, ‘ 


2 h 
h+2X 


n-1 

^ K(n-v,0) 0(vh) 


A calculation carried out by hand using X g /h**3 and X s *.05 shows close agree¬ 
ment with previously calculated approximate solution, 'Hie equation In the form above 
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has been programmed for Whirlwind. The first results which agree with the hand calculated 
curve were recently obtained. 


Joseph B. Walsh 


COULOMB WAVE FUNCTIONS 


The programs for both the regular and irregular coulomb wave functions for L*0 
are now complete and an output routine is now being written. Since this is the last 
formal progress report for the Project, we will summarize in some detail the nature of 
the computations involved, reviewing some of the work previously reported. 

RANGE OF PARAMETERS 
The ranges of parameters and L will be 


in the following sense 1 . Knowing the value of T{?jp) and F* at P. , we 

values of F and F' at p + by using the Taylor expansion J* 

For machine computation, It la more convenient to uee the notation 

> ia r> a 

riP+Ae). £ C7 Xn) ; r< f+Jf) - L- 2 ncr<n) 

'• ; n-O J' ' A] n=0 


find the 


0 4 > i 12.5 


5 r 4 


with intervals 


Af *A>> » 0.2 

A? might possibly bo changed to 0.1 but thle Involves only a trivial modification of 
the present program. 


METHODS OF COMPUTATION 


a) Regular Function 


For the interval 0 £ ^ £. 1, the power series is used 

VAr> - c .<>> r £ B n 

/ n=0 

C o < > ) " c?^-! 


B o * B 1 


(irf-1) (irt“2)B 




About twonty terms are needed to give eight digit accuracy for the range of P and *> 
given above. In the Interval 1 < ft 25, .e us. numerical Integration of the dlffe 


equation 


differential 


.. 1-29 

F" + (‘ p J . 0 


gXl) = F-( f ) 

cr* 2) - - |(i- F( j>. ) (up 2 

f> n(nfl) <T <nfl> + (n 2 -n) (T (n) + if- J 2 0 in ' 1) 


+ ^ 3 <r (n ' 2> = Oi n^2. 


Enough terms in (3) are taken so the number of significant figures in tbe new F and F' 
is the same as that in the old F and P’. Of course, near a zero of F or P', merely 
taking a lot of terms is not sufficient to preserve the number of significant digits. 
To show this we write, 




A « 1 - ; ( 




dp* a - y + 


; B = ( J/”) + ... 


where A and B are infinite series determined by (3). Lot us assume that A and B a-e 
known exactly and that F and F’ are related to the true values of the functions, F ( and 
F t ’ , by 

F - F t (1 + 

F 1 = F t ’ (1 +< 2 ) 

lej - IV - € 

We have therefore 
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r. S ) « F t ( +Af> >) + f/A r t ( ye» > ± Bf t ' (JJ )/ (5) 

where F c star.is fc- •- calculated value of P. Since A F t and B F*' are of the order of 
unity , it if possible for the error tern in (5) to bo of the order of £ . But if 
Ft ) ll 1 r. magnitude than AP( P ) or BP'( fi ) because of 

strong subtract: we will evidently lose sane significant figures. In the computations 

performed so far tie strongest zeros encountered were of the order of 5(10)” 3 . At those 
points, we aigfct expect only five significant figures if seven are carried at other 
points, as is the case ir. the present computations. The machine tine for calculating F 
an d P 0 ' for a single * ? varies between 400 and 800 nilliseconds depending on the number ° 
of terms used. For most of the range in f> (2 t f* i 25) it is sufficient to take terms 
up to n = Q in (3) with a corresponding machine tins of about 500 milliseconds. 

b) Irregular Function 

"Hie irregular functions G 0 and G 0 ’ are calculated by the technique indicated 
in (3). The starting points for the procedure are the transition line p = 2 £ for 
) * and the points j* - 4 f c' 6 *7 * 2^ for the remaining region. For f* = 2 {2 £*?£{} 

and p m 4 t o 4} • 2, we use the integral representation for G c and G 0 ' 2 . 


F - -K).91794506 I 

° ( 

F * - 40.33103212 / 
o 

b) Irregular Function 

f = 20 
> - 5 

G q * 4*1.1657166 Numerical integration from 

p « 5, jO «: 10 

G s 4-1.1657160 Asymptotic formula 

o 

References 

1. Abramowitz and Rablnovltz, Phys. Rev. 96, 77 (1954). 

2. Newton, Atomic Energy of Canada, Limited, CRT-526 (1952). 


Integral representation 
or steepest descent 


G o<£>> * 

G o= A c < > ) //' + (6 > 

VV > - H i~ 

*o - s - ip( ^ 2 y 


Gaussian quadrature is used to evaluate the definite integrals with the cut-off in ? 
chosen so as to give seven to eight significant figures. ' 






J 7 ' we Ea J' u8e the method of steepest descent 1 

integrals and obtain a rapidly converging serlos in Inverse powers of>.‘ The method 
of computation is subject to the same limitation discussed in reference' to the regular 
f action, namely that near the zeros of the functions the number of significant digits 
accurately retained is (7 - 2) where Z is the number of zeros between the decimal poin 
and the first significant digit. p 



SAMPLE COMPUTATIONS 

Regular Function 


/ = 10 


> ‘ 5 


F = +0.91794503 ) 

° r 

numerical integration 

F - +O.J3103219 1 

o / 

fro " p “if" /• 


Arnold Tubls 
Aaron Tomkin 
Zoltan Fried 


FIRST APPROXIMATION SOLUTION ON ORE BODY 

Proceeding from the favorable results reported in the previous Quarterly 
Progress Report No. 19, the investigation of the first approximation solution was di¬ 
rected towards the determination of the optimum size and shape of the ore body to be used 
in the inverse problem of Interpretation of field data. However, certain difficulties 
were met which have not been completely solved, specifically: the physically modeled re¬ 
sults have been found to be somewhat unreliable, and a paradox has been met in the deter¬ 
mination of the size and shape of the fundamental building block (optimum ore body) for 
the inverse problem. The character of the modeled results is always consistent but the 
numerical values are not, and hence more accurate experimental work will have to be done 
before a final decision can be made concerning the feasibility of the first approximation 
solution. Also the paradox will have to be resolved as regards the optimum ore body size 
and shape. 


The first approximation solution does have its merits and as such will be of 
valuable assistance in the interpretation of field data and the economic evaluation of a 
possible mineralized area. 

In conjunction with course 6.538 (Electronic Computational Laboratory) a com¬ 
plete theoretical and numerical solution to the three-vertical-layered earth has been 
done utilizing the IBM Type 650 Magnetic Drum Computer. That is to say, the apparent re¬ 
sistivity profiles for a vertically layered earth have been obtained for both pole and 
dipole senders for a large number of geometries and resistivity contrasts. The programs 
for these calculations are available for further work, and the results obtained are pre¬ 
sented in the report VERTICAL LAYERING IN GEOPHYSICAL PROSPECTING submitted as a term 
project in course 6.538. The rosults obtained have been of assistance in the evaluation 
of the first approximation solution and also as a tabulation of resistivity profiles for 
certain idealized geological situations. These can bo used to assist in the Interpreta¬ 
tion of field results before further investigation on hijh-speed computers is made. 

Norman F. Ness 
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DYNAMIC RESPONSE OP SHEAR WALLS 

A "shear wall carries loads applied In the plane of the wall parallel to the 
length of the wall. Walls In actual buildings carry portions of transverse loads (wind, 
blast, earthquake, etc.) In this manner. 

In this investigation the actual wall is replaced by a dynamic model consisting 
of a number of concentrated masses Inter-connected by weightless springs (see figure). 

The properties of the springs are obtained by equating the deformations of the actual 
wnll and dynamic model under similar stress conditions. With the spring properties 
known, the equations of motion are obtained for both horizontal and vertical motion of 
each mass point. 


RESPONSE OF A PIVK-STORY PRAME BUILDING TO DYNAMIC LOADING 

'Hie object la to determine the response of actual multistory rectangular 
rigid-frame buildings subjected to blast loads. A more complete description may be 
found in the March 1956 Progress Report. 

The control program Is written and being debugged". The load program is 
written and debugged, as is the Integration program. The resistance program is In two 
parts. The first part, the solution of a matrix of the slope-deflection equations, Is 
currently being tested In conjunction with the second part, which reviews the matrix 
solution for changes from elastic-to-plastlc or vice-versa and modifies the matrix ac¬ 
cordingly. This last part has been successfully tested by itself. On completion of this 
test, the various programs will be combined and tested for complete solutions. 


Throe different size walls are being studied for two dlfforont load conditions 
and two different sized grids for a total of twelve solutions. The solutions are re¬ 
stricted to the elustic range. The number of equations varies with the wall size and 
fineness of grid; the minimum number is 15 for a wall 4’-0" square with a 2'-0” square 
grid, and the maximum number Is 77 for a wall 4*-0" by 12' -O" with a l'-4" grid. The 
equations of motion are solved by numerical integration. Equations of the form 


x 


n+1 


* 2x 


n 


x , + X (At) 2 
n-1 n 


whore X |l+1 displacement at n+1 st time interval 
x^ = velocity at nth time interval 


Ralph G. Gray 


CALCULATIONS FOR SPHEROIDAL NUCLEI 

Tlie computations which were originally projected under these problem numbers 
is now complete. Before submitting a final report, however, two additional investiga¬ 
tions will be undertaken. 

a) Problem 235, Bound States in a Spheroidal Well. It would bo desirable to 
try to include the effect of a spin-orbit coupling term upon the energy eigenvalues. 

Such a calculation will be carried out if it doos not involve any substantial revision of 
the existing programs. 


and 


x =■ acceleration at nth time interval 
At & length of time interval 


arc used to predict deflection at time . A second equation of the form 


X (At) J 


x n+l * *n + i n <At) + "hr 


W A *> 


b) Problen 319, Scattering from a Spheroidal Square-Well of Zero Energy 
Neutrons. There Is a relatively simple way of approximating the case whore the scatter¬ 
ing potential is complex. In essence it Involves the calculation of the Integral of the 
square of the wave function (for a real potential) over the volume of the scatterer. A 
program to do this has been written. 


Jack L. Uretaky 


FINITE BENDING OF THIN, SHALLOW SPHERICAL SHELLS 


is used to correct the previous estimate. An iterative procedure is required to obtain 
a stable solution since x is Involved. After a set of consistent displacements at 
time t has been obtained, the principal tensile stresses at a discrete number of 
points on the actual wall are calculated. Those are compared with the cracking stress 
of reinforced concrete. If cracking has not yet occurred, the process Is repeated. At 
the end of each time interval, horizontal and vertical displacements, stresses and re¬ 
actions are displayed on the scope. 

The proper value of A t la obtained by comparing solutions using several differ¬ 
ent values of fit. The program is currently In production. 




Computations of the bending of a uniformly loaded shallow spherical shell have 
been completed in the case of simply supported edge. In particular, buckling loads for 
very shallow shells have been obtained. The results show that the convergence of the 
employed power series method as outlined in the last Progress Report is sufficiently 
rapid as to give fairly accurate solutions for stable and unstable equilibrium positions 
in tho post buckling region, provided the shell Is sufficiently shallow. Calculations 
for less shallow shells indicate, on the other hand, that solutions can be obtained in 
this framework only for loads below the buckling load and in some cases also approximate 
solutions beyond buckling. The power series method seems to be able to reveal equili¬ 
brium states close to the buckling region in the latter case. Similar phenomena were 
observed in the case of clamped edge. Two circumstances may be mentioned here that are 
responsible for the break-down of the method. In the region of buckling the vertical 
deflections of the shell show a complicated pattern which cannot be approximated by a 
polynomial of sufficiently low degree. The second resson is that the initial estimate 
is not sufficiently close to the actual solution; hence the number of iterations to be 
carried out to obtain a solution would be extremely high. 

The results obtained during this quarter are now being evaluated and analyzed In 
somewhat more detail with respect to some special solutions that have been obtained 
hitherto for this problem. A more complete description of th# results will be given in 
the next Progress Report 
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THE NUMBER DISTRIBUTION Of ELECTRONS AND PHOTONS 

The Whirlwind coding of the grid doecrlbed in Progress Report No. 13, p. 55, 
has been successfully carried out. At prosent, » program selecting paths according to 
the Monte Carlo method described in Progress Report No. 18, p. 12, is being tested in a 
simple case (n*m=5), and will be easily extended for larger n and m 

The next step which does not Involve any logical programming will be the com¬ 
putation of the products of transition probabilities and convolutions of exponentials in¬ 
volved in the formula given in Progress Report No. 14, p. 16. It will be partly carried 
out by Dr. B. Rankin, and continued by myself starting on August 1st. 

Raymond F. Stora 

ASYMPTOTIC SOLUTION OF A DIFFERENTIAL EQUATION 

Research lias been continued on the solution of the differential equation with 
turning points’ described in the last report. The third order equation 

u'" + A 2 Z u 1 + 3(iA 2 u = 0 

X large 

p constant, complex 


which is a special case of the previous one, is now being solved by use of the Laplace 
Integral. In particular, the asymptotic expansion of the solution 



c(Ai) 

A comparison will then be made with the asymptotic series expansion obtained by 
Rudolph E. Longer with solutions computed numerically. 


Joseph Rershenov 

AN APPLICATION OF MONTE CARLO METHODS TO NEUTRON DIFFUSION 

airing the past quarter all routines have been rewritten and retested. Several 
errors were found. The first assembly of the individual routines has been run success¬ 
fully. Tests of problems, for which the analytical solution is known, are continued 
This problem is programmed and coded by M. Troost, Research Assistant in Chemical Engi¬ 
neering. 

Marius Troost 
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CONDENSATION IN TUBES SUBJECT TO THE EFFECTS OF VARIABLE VAPOR VELOCITY 

Two runs at different conditions of entrance vapor flow but with the some vapor 
and condensate properties (those of water at 212°F) appear to give reasonable results. 
However, there has not yet been time to check these results with existing data. 

The stopping criterion in these cases was the transition of the condensate from 
laminar to turbulent flow. The program for the turbulent vapor - turbulent condensate 
caBe is being completed in order to continue the laminar condensate runs to the condition 
where all the vapor in the tube has condensed. 

It appears that the turbulent vapor-turbulent condensate case analysis cannot 
be simplified to contain only throe dimensionless parameters describing the geometry, 
initial conditions and vapor and condensate properties as was possible with the turbulent 
vapor-lamlnar condensate case. 

Jukka A. Lehtinen 

RESPONSE OF A SINGLE STORY REINFORCED CONCRETE BUILDING TO DYNAMIC BLAST LOADING 

The object is to study the dynamic response of the roof of a single story 
reinforced concrote building for various loading functions. The roof is a system of In¬ 
finite degrees of freedom consisting of a one-way slab supported on relatively stiff 
beams. This is resolvod into a two-degree system by assuming deflection shapes for the 
one-way slab and its supporting beams and using LaGrange's equation to obtain two dynamic 
equations for acceleration of the beam and slab. These are Integrated for deflections by 
an "acceleration impulse method" which utilizes backward differences. 

The program is currently in production. 

Ralph G. Gray 
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Part II 

Project Whirlwind 
I. REVIEW AMD PROBLEM 1MDEX 


This report covers the specific period of March 20, 1956 to June 10, 1956. 
IXirlng this time, 95 problems made use of 398.03 hours of the 514.25 hours of Whirlwind 
computer time allocated to the Scientific and Engineering Computation (S&KC) Group. 

Of the 514.25 allocated hours of computer time, .2% was down time due to computer 
malfunctions. The remaining 104.62 hours of the allocated time were used for terminal 
equipment testing and calibration, demonstrations, and various inter-run operations not 
logged to specific problems. 


rhe 95 problems run during this quarter cover some 20 fields of applications. 
The results of 32 of the problems have been or will be Included in academic theses. 

In these 32 problems, there are represented 25 Doctoral theses, 3 Naval Engineer's, 

5 Master's, and 4 Bachelor's. Twenty-six of the problems have originated from research 
projects sponsored at MIT by the Office of Naval Research. 


Two tables are provided as an index to the problems for which progress reports 
have been submitted. In the first table, the problems are arranged according to the 
field of application, and the source and amount of time used on WW1 are given. In Table 
2-II, the problems are listed according to the principal mathematical problem Involved 
In each. In each table, the letter after the problem number indicates whether the 
problem is for academic credit and whether it is sponsored. The code is explained in 
section 2.1, Introduction. 
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PROBLEM INDEX 


PROBLEM 


Field 

Aeronaut ic.a I 


Aerupitia lea 
>trivarcu Uruup 


Ueacrlj»l>uu 

Transient taa|i«r*tur« uf a Im>i t ,p*< !»*• 

liana i an I r««|><aiaa u/ aircraft structure# to aerodynamic *•■«' 'W 
Poet fa Ilure reapunav of aircraft structure# subjected tn blast 
Ida itlnc 

I'araaatric atuil> uf coupling ami damping 
Vartlcal tall Itteda <iua to rolling pull-up 
Itlast response of rotor blades 

tracuai of stability derivatives turn flight leal da la 


Miami cal 
Aiik Inaar Inn 


rranauni<a la dlsilllaitun column# 

Critical aaia calculations for crl Intfrlctl (MMln 


Civil 


Kna JD«*ar Inc 


Ueap,ms«> of a single aton coocral* building 
U“al«n uf spherical shell segments 


Nlectrleal 
Engineerlot 


a«iar*> 

Conversion lab. 


Linear progroaming 

Production lor transputlatum problem 
duatlulluo orrur 
Propagat ion uf ruund-off arror 
Uptlalutliifl uf a 1 1 area tor .-nntrol »it«a 
Opt la lost loti of raa-atr cooling ayataai 


U*olo|) founrr synthesis for crystal atru<stur«a 

Inatrumen tat ion Data reduction for X-l fire control 
Laboratory 


Lincoln Labor*t><rv Eigenvalue prcMiloai for propagatlun of alact ruaagnet lr a»»«a 
General raydlat solution 
Tropospheric propagat loti 
frror aaalyaia 
Prediction analysis 

Atmospheric propagat too of radio wave* 

Coverage analysis 


Mechanical 

Engine*ring 


Meteorology 


Naval 

Architecture 


Thermodynamic and dyntaic efforts of eater Injection Into high 
temperature high velocity gas atreaaa 
Laetnar boundary layer of a a toady, *-<*preaat bl« rl«* la the 
entrance region of a tut** 

Hulling bear Inga 

Plater -tablllast ton on a heated plate 
lleaa vibration 

Trane lent hoat flow tn eoltda 
Matrla Multiplication 

Aayaptotlc Integration of equations concerning toroidal shell 
Investigation of the vorttelly field in the general clrculatloo 
of the ataoaphere 

Hpo, iral analyala of ataoaplierlc data 
Statist leal and dynaaic aethoda In foracaattng 
Weather prediction 

Computation of variances and covariances 
Solution of transverse eeb t rise 


Naval Supersonic 
Laboratory 
Nuclear 
engineering 


Office of 

Statistical 

Services 

('!>) a Ira 


te rv rase chan 1 was 
Laboratory 


Diffusion boundary layer 

Calculations for tho MIT reactor 
Non-uniform fuel distribution 
Uatormliuitlon of critical taaaa 
flui leveling in homogeneous reactor - Part 1 

finding eigenvalues of an uyaMirtc Matrix 
Temperature distribution In a bean 

Determination of pliaae ahlfla froa aXpertmenlal rroaa-eocttuna 

Atoailr integrala 

Theory of neutron read Inna 

Scattering froe oxygen 

'll’* *a applied to face and body-centered iron 
Energy levels of diatonic hydrtdeu 
•(valuation of two-center Molecular Integrala 
Analysis of air atiueer data 
Knergv levels la dlatoailc hydrides L1H 
Atoailc mu vo functions 

1‘olarlaabll lty effect a In a time and molecule* 

Puro and tnpure potassium chloride crystal 
imta reduction 


Spectroscopy Ub Complex spectrua analysis 


"roMem Mliiulea t 
•<uabei VH I 1 tee 

170 C. 03.3 

j>. c. 13.a 

330 C. 170.a 

334 C. Ml.3 

358 H. 33.3 

364 C. 00.S 

317 C. 262.0 


241 S. 100.0 

270 H. 996.0 

.134 0. 110.4 

373 B. 3.0 

2 ia 337.1 

336 C. 60.3 

335 0. 30.0 

337 B. SO.O 

364 C. 150.3 

336 C. 139.7 

261 C. 296.2 

244 C. 6S.3 


103 L. 547.3 
272 L. 173.• 
300 L. 220.2 
312 L. 237.1 
337 L. 306.7 
371 L. 237.4 
377 I.. .70.4 


120 li.N. 167.7 

190 N. 60.S 

293 C. HS.o 

333 A. 74.4 

330 A. 66.H 

342 B. 1M0.3 

343 R. 70.0 

J63 A. 14.3 


Supervisor or 

J. C.Loris 
L.C.Schalt 

H.U'laite 

K. Vetmore 

II. Parechantan 

K. foe a 

L. L.Maaxola 


S.K.Uavta, Jr. 

J. R.Poeell 
B. Landry 

K. Traum 
|i. Arden 
J.Dennis 
II. Pop Ip 
A.I.Green 
J.Dennis 
tt.Moronoy 
a.J.Buerger 
J.Stark 


II. B. Height 
0.< .Sponaler 

H. B.Height 
L.i’oteraon 
L. Peterson 
• .Mason 

I. . Potereon 


A.Erickson 

T.V.Toooi 
A.Shaabaty 
T.Y.Toong 
8.H.Crandall 
H.Gelman 

R. R.Archer 

S. II.Crandall 


226 0, 76.g 
306 0. 217.8 
341 C. 604.0 
343 C. 633.4 
330 D. 203.9 


D.Cooley 
R.Salliman 
•-Holley 
J.Austin 
D.Oilman 


330 B. 85.4 C.Brandt 


307 B. S3.2 J.Baron 


M6 A. 73.3 T.Cantwell 

351 B. 126.7 A.Sutter 

367 B. 123.3 J.Barnett 

373 B. 17.0 ft.Kennedy 


•W5. 134.5 r.M.Verruh 

3«0. 2BH.S f.M.Verxuh 


162 H. 36.4 

234 N. 2.3 

245 N. 2059.0 

246 H N 279.3 

333 N 127.1 

280 H. 349.1 

262 N. 700.1 

272 N. 243.9 

278 N. 673.1 

*88 N. 1670.8 

290 N. 667.6 

309 U,N. 276.9 


F .KppWng 
H.K.Neabet 
H.Pesbbech 
A.Twmktn 
J.Wood 
A. P reman 
A.Agha, uilai. 
P.Schc .i 
A.knro 
H.K.Neabet 
• 

L.l'. Howl and 


136 C. r.7ft.l D.T.Koaa 


346 B. 123.3 J.Lindner 




MntbMat teal >rob I am 


1. Matrix algebra and equal Iona 

(tout uf determinants! equal lull 
Linear equal tuna 
Slgenvnluea 
Bigeuvaluea 

Orthonormal I tat tun, elgenvaluee 

Inversion 

Linear equatluiis 

Kquat tuna 

Ordinary algobra 

Diagonal last Ion, Inversion 

Inversion 

Arithmetic operatluns 

Solution of 5 «leu 1 1 annum firal degree aquations 
Inversion, eigenvalues 
Eigenvalues, (act ,Mug polynomials 
Arithmetic operation* 

2. Ordinary differential equations 

Sewn mm-linear first order 
System 

Non-linear firal order 
Pave equation 
Non- I Inear 

Set of differential equations 
Slmul tanmaia dlffeianttal pquatluna 
Second enter 

Nor linear differential equations 
Second order 

Three atm.il raneuua aecond order 
Pair of simultaneous aecund order 
Ordinary differential equations 
Set of differential equal Iona 
Second order 

J Partial differential equations 
Second order parabolic 
Plrat order system 
Second order parabolic 
Fourth order 

Second order, variable coefficients 

4. Integration 

Integral tranalunuttut, 

Kvalnation of integrala 
Overlap Integrala 
Ovsrlap Integrals 
Integral lone 
Fmane I Integral 

Integral Iona 
Integral Inna 
IntegralIon 
Integration 
Integral Inna 

3. Transcendental equal lone 

Curvn fitting 
Non-liuear equations 
Non-linear equations 
Curve fitting 
Non-linear equation* 

it. Cumpiex algebra 

Complex roola and function evaltottion 

7. tt»ta reduction 

Data reduction 
Mia* dJetrtli.it ion 

8. Fourier aertea 

Fourier aynthoaie 
Fourier ayntheel* 

9. Linear programming 

Lieoar programming 

10 . Itetietlce 

Nvmluatian of covariance and variance 

Autocorrelation 

Mean and second womeni 


UeraiiuH 266 A. 

Croat 370 B. 

Ml agonallasI ton 27S a, 

diagonal Iasi ion vaa N. 

Scheldt process. dsaguaalI patina 3B" *• 

(Ihkii 312 l. 

Croul 317 C. 

Sairla iteration 330 C. 

Iteration 338 C. 

•ppmatmatinns, Croul 341 C. 

rout •< • 

in rev 4 346 B. 

Crout 339 B. 

Croat, iteration 344 C, 

Frane. Hitrhrock 363. 

Direct .177 L. 

Fourth order Kutra-GIII 120 U,N. 

i i l i IB 9 N. 

Sec (aid .,rder Kutta-Cill 241 b. 

Milne predlctort.,rrector formula 243 N, 

0111 297 II. 

ant 3i.1 A. 

Finite difference 334 C. 

Power eerie* 342 B. 

Pntturbutton procedure 34.*, ft. 

Finite difference 331 B, 

Traps got da l rule 1M B. 

Finite difference Mil A. 

Finite difference 3*4 C. 

Kittle Oil l 371 L. 

Pintle difference 374 B. 

Finite difference 179 C. 

Finite difference 22* 0. 

Finite difference 236 C. 

Finite difference 13# A. 

P'nlte dtfferencu, trepexoldal 

integration 37;i B. 

Algebraic recursion formula 234 N. 

Approxlestion by algebraic functions 246 B,N 
Rvalue11on of analytic force 260 N. 

(valuation nf analytic forme 262 N. 

Slapeon** rule 378 N. 

Converalnn p>>e« r ear lea with a 

complex argument 300 L. 

Barnet and Coulaon expansion 300 B. 

Uauaa quadratura 312 L. 

rihite di rference 

I'rapecoldal rule 369 B. 

Series aubattcutloa 372 B. 

beset squares 163 N. 

Steepest descent 264 C. 

Iteration 279 N. 

Least squama 373 N. 

Steepen! descent 3V3 C. 

Iteration IBS L, 

Polynomial fitting, etc. 136 C. 

Arithmetic operatligis 244 C. 

Direct araluetlon 361 C 

Direct evaluation V» 0. 

Staple* relaxation 319. 

Direct 160 D. 

Direct 166 S. 

Diracl 357 II. 


Table 2-1 


Current Problems Arrange.I According to Field of 


Application 


Table 2-II Current Problems Arranged According to the Mathematics Involved 
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2. WHIRLWIND CODING AND APPLICATIONS 


100. COMPREHENSIVE SYSTEM OF SERVICE ROUTINES 


2,1 Introduction 

Progress reports as submitted by the various prograwirs are presented In 
numerical order In Section 2.2. Since this sumaary report presents the combined efforts 
of DIC Projects 6345 and 6913, reports on problems undertaken by sesbers of the Machine 
Methods of Casiputatlon (MMC) Group have been emitted fro* Section 2.2 of Part II to 
avoid duplication of Part I. For reference purposes, a list of the MMC Group problems 
appears below. 


Letters have been added to the problem numbers to indicate whether the 
problem is for academic credit and whether It Is sponsored. The letters have the 
following significance. 

A Implies the problem Is NOT for academic credit, Is UNsponsored. 

B Implies the problem IS for academic credit, Is UNsponsored. 

C implies the problem Is NOT for academic credit, IS sponsored. 

D Implies the problem IS for academic credit, IS sponsored. 

N Implies the problem Is sponsored by the Office of Naval Research. 

L Implies the problem Is sponsored by Lincoln Laboratory. 

The absence of a letter indicates that the problem originated within the SAKC Group. 


List of Machine Methods of Computation Group Reports 


The following problems used Whirlwind I computer time during this quarter, 
but are reported on In Section 2.2 of Part I. 


122 N. 

172 B,N. 
203 C. 

231 B,N. 
240 A. 
319 B,N. 
329 N, 
337 N. 

361 B,N. 
368 B,N. 


Coulomb Wave Functions 


A. Temkin 
A. TubIs 
2. Fried 


Energy Bonds in Graphite 

Rosponse of a Multi-Story Frame Building Under Dynamic 
Loading 

Reactor Runaway Prevention 

Number Distribution of Electrons and Photons In Cascade 
Zero Energy Scattering Cross-Section of a Spheroidal Well 
First Approximation Solution on Ore Body 

Non-linear Second Order Differential Equations In the Theory 
of Elastic Shells 

Growth of Fatigue Cracks 

Condensation In a Vertical Tube 


F.J.Corbato 

R.G.Gray 

M. Troost 
B. Rankin 
J. Uretsky 

N. F.Ness 

H, Welnltschke 
J.B. Walsh 
J, Lehtinen 


Several minor and obscure errors have been detected In the C8 conversion 
program. For the most Important of these, the necessary corrections already have been 
Incorporated In the program. This corrects the situation In which the current address 
Indicator could be set definite" even though It was set by a word containing the syllable 
r when r Is indefinite (see C.8.Manual, Chapter XIV). The conversion post-mortem 
has been printing out the wrong register as the address for Program too long at " when 
a ditto block runs beyond the beginning of the PA. The necessary changes have been dis¬ 
covered and checked out, but due to the difficulty of re-recordlng magnetic tape unit 0, 
the correction has not yet been made permanently. A logical error In the program which 
gives the wrong result only In a peculiar and quasl-legal situation, a word like ts5r-" 
after an "1.NR ' In which the minus sign Is taken as preceedlng the "r", has been found 
and a correction has been checked out. 

The manusl modes of the new Generalized Post-Mortem (see Summary Report *45) 
are being checked out. Most of the errors hsve been discovered and corrected. The new 
program will display memory on the scope at a rate comparable to that of the scope post¬ 
mortem In the Utility Control Program, a considerable advantage In speed over the present 
program. Besides, not only zeros, but any block of Identical lines Is suppressed and 
only the first and last llnBs are given 

The fp tape reading program and the fled nodes sill be checked out after the 
manual modes are operating satisfactorily. As the Utility Control Program must be 
modified to accommodate the new Generalized Post-Mortem and some other changes seem 
desirable to take advantage of new in-out equipment, etc., a largo number of revisions 
seem necessary. The new system will be working in the near future. 

M. R. Weinstein 

Digital Computer Laboratory 


120 B , N , THE AEROTHERMOPRESSOR 

During the past quarter, computations were carried out to check wltb the 
experimental results obtained from the Aerothermopressor In the Gas Turbine Laboratory, 

The agreement Is very good. More experiments and theoretical computations are to be 
carried out In the next quarter. 

A. Erickson 
Mechanical Engineering 

126 C. DATA REDUCTION 

Problem 126 IS a very large data reduction program for use In the Servomechanisms 
Laboratory. The overall problem Is composed of many component sections which have been 
developed separately and are now being combined Into complete prototype programs. Des¬ 
criptions of the vsrious component sections have appeared In past quarterly reports. After 
the development and testing of the prototype Whirlwind programs Is completed, the programs 
will be re-coded for other, commercially available, large scale computers, (probably the 
ERA 1103, IBM 701 and IBM 704 computers), for use by Interested agencies for actual data 
reduction at other locations. The programs are currently being developed by Douglas T. 
Ross, ttavld F. McAvlnn, Walter E. Weissblum, Benson H. Relief f, and Dorothy A. Hamilton, 
Servomechanisms Laboratory staff members with the assistance of John F. Walsh. This work 
Is sponsored by the Air Force Weapons Guidance Laboratory through DIC Project 7138. 
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The nature of the problem require* extreme autonaticlty and efficiency In the 
actual running of the program, but alao require* the preaence of human operator* In the 
computation loop for the purpoae of declalon making and program modification. For this 
reason extensive use Is made of output oscilloscopes so that the computer can communicate 
with the human, and manual Intervention registers so that the human can communicate with 
the computer In terms of broad ideas, while the computer Is running, and have the computer 
program translate these Ideas into the detailed steps necessary for program modification 
to conform to the human operator's decision. The program which does this translation 
and modification Is called the Manual Intervention Program (MIV). The most recent version 
of the prototype data reduction program is called the Basic Evaluation Program. 

During the past quarter work has continued on the Evaluation Program and the 
MIV Systems. A number of changes and revisions wore made and the complete system has 
been used to make several test runs for comparison against previously computed values. 

The elaborate scope plot which was described in the past quarterly report has been written 
and incorporated into the manual Intervention system so that all pertinent functions 
within the evaluation program may be plotted for any range specified by throwing appropriate 
switches. 

As mentioned in previous quarterly reports, the programs which are being devel¬ 
oped on the Whirlwind Computer are prototype programs which are to serve as models for 
production Univac-^cientIfic 1103 computer programs. Several meetings have been held 
with personnel of the Air Porce Armament Center, Eg!in Air Porce Base, Florida, for the 
purpose of coordinating these efforts with their 1103 installation. Besides the design 
of these programs, procedures have been Initiated for the procurement of additional Input- 
output equipment for the 1103 computer for the use of this and other problems. In particular, 
a Charactron display tube for visual and photographic facilities, and a complete manual 
intervention system have been ordered for the 1103. 

A proposal has been made and accepted for the addition of a Flexowrlter key¬ 
board Input facility to the Whirlwind Computer as a result of the requirements of problem 
126. This facility will prove highly useful to other users of the Whirlwind Computer as 
well. The modification, which is scheduled for completion In July, consists of assigning 
new si addressee to the Flexowrlter keyboard In test control (and also an auxiliary Flexo- 
writer at the manual Intervention console used by problem 126) so that the Flexowrlter 
keyboard may be used in the same way that a photoelectric or mechanical tape reader now 
operates. In this mode arbitrary alpha-numeric information may be typed directly into 
the computer through the ln-out system, and arbitrary interpretation of this information 
programmed Into the computer. The Comprehensive System will be modified so that this 
facility may be used for ordinary CS input when the occasion warrants. Due to the low 
time efficiency and reliability of this mode of input it Is expected that its use will 
be restricted to applications of a manual Intervention nature. Applications of interest 
to the general user, however, will be to operate the Director Tape Program for program 
manipulations and operating demonstration programs, of major Interest to problem 126 
is the fact that individual keys on the typewriter may be treated as separated intervention 
switches so that this new facility in effect will add a large number of intervention 
switches to the present system, in addition to the more general flexibility mentioned 
above. 


In the past quarter, a brie* study was made of the poaalbllity of Incorporating 
a new Instruction Into the Whirlwind repertoire to allow a type of microprogransalng on 
the Whirlwind Computer. The proposed l.iat.uctlon, ml. would allow referencing of lndlvldusl 
command pulee output (CPO) lines eo that special purpoae lnetructlons could be manufactured 
by the program*. The study was restricted to consider schemes which would require that 


a minimum of additional hardware be added to the Whirlwind Computer and that the normal 
operation of the computer would be unaffected. These conditions required that the new 
facility be built out of existing CPO units and consist primarily of a rewiring of certain 
parts of the control element. Since the program register, PR, is still a part of the 
hardware of the computer although it is not used at present, this register was chosen as 
the focal point for the ml instruction and an assignment of CPO units to the 48 gate tubes 
associated with the PR circuitry wus made. It was found that a simple, logical circuit 
could be made to load the program register from arbitrary core memory locations and use 
the word thus obtained In PR to energize gate tubes to control individual CPO units on 
selected time pulses. It is also possible to have an arbitrary number of addresses In 
the ml instruction, the meanings of the addresses being controlled by the selection of 
the proper CPO units. The design which evolved calls for the execution of one ml 
instruction to automatically Initiate the reading of an arbitrarily long sequence of 
core memory address locations Into PR for execution one at a time. In this way, arbit¬ 
rarily long and complex instructions can be manufactured. 

After the design was completed, however, it was found that the flexibility which 
appears to result from the system is in fact quite illusory due to the very tight design 
of the Whirlwind logical-arithmetic element. A major drawback seems to be that the only 
access to the accumulator in the Whirlwind arithmetic element is through the A register, 
so that the AR is always tied Inflexibly to the AC. The net affect of this is to require 
an inordinately large amount of extraneous shuffling of information in order to use the 
mi facility so that the system is extremely inefficient. Although one of the basic 
premises upon which this study was undertaken was that the system would not be required 
to be efficient time-wise, it turns out that it is almost always easier and faster to 
write a Whirlwind program in the existing order code, to do a particular Job, than to 
manufacture an ml instruction sequence to do it. For this reason the present design will 
not be given consideration for an actual proposal. 

At this time it Is very difficult to estimate how much this picture could be 
Improved by allowing the addition of several new and more general CPO units, but it seems 
probable that quite a drastic revision of this evaluation of an ml facility could result 
from the choice of a few appropriate additional CPO units. The details of the present 
study will appear in a Servomechanisms Laboratory memorandum, and a supplementary study 
to consider alternatives may be made In the future. It is felt quite strongly that not 
enough Information is known at this time about what a microprogramming facility should 
consist of, nor how one would be used If available, so that some inexpensive facility 
should be designed and used before consideration Is given to building an actual micro¬ 
program computer. It may be, however, that such an experimental facility can best be 
achieved by simulation on existing computers rather than rewiring of existing computers. 

D. T. Ross 

Servomechanisms Laboratory 


162 N. NUCLEAR SCATTERING PHASE SHIFTS 

This quarter has seen further work in fitting the phase shift combinations to 
the experimental curve at the first minimum. It has been necessary to take a much finer 
mesh in this region because of the sensitivity of the phase shift eqjatlon. 

Future plane Include further work along the same lines to fit the curve beyond 
the first minimum in such a way that the results can be projected to encompass the 
values of the entire curve at high energies. 

Programmers working on this problem are K. Mack and E. Campbell. 

E. Mack 

Nuclear Science Laboratory 
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168 D. INDICIAL DOWNWARH BEHIND A TWO-HIMKN8IONAL WING 

ThlB problea is now completed. The result* will be found In an Sc. D. thesis, 
entitled Indlcial Downwash and Its Structural Effect on the Horizontal Tail , submitted 
by N. P. Hobbs to the MIT Aeronautical Engineering Department, June 1956. The results 
of the problea will also be published as a Wright Air Development Center Technical Report 
56-164. 

N. P. Hobbs 

Aeroelastlc and Structures Research Lab. 


known to be impractical or abaurd. Also, provision has been made for storing the program 
and basis table on the buffer drum, thus making the whole auxiliary drum available for 
storage of cost data. With these improvements, the program has the capacity to handle 
up to approximately 10,000 significant costs. 

In the course of trying to solve a large problem (m = 60, n * 291), a number 
of coding errors showed up that had not been detected earlier. These have been tracked 
down and the 60 x 291 problem has been solved in 30 minutes of computer time. 

J. B. Dennis 
Electrical Engineering 


179 C. TRANSIENT TEMPERATURES AND STRESSES IN A BOX-TYPE BEAM 

The programs have been modified to obtain transient temperature and stress 
responses for the same structure for different thermal Inputs. 

During the next quarter, the program will be repeated for additional thermal 
inputs to define the dependence of these responses in terms of input parameters. 

J. C. Lorln 

Aeronautlcsl Engineering 


226 D. CIRCULATION OF THE ATMOSPHERE 


IXirlng this quarter a program has been written by Martin Jacobs for solving 
by relaxation the equation 
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193 L. EIGENVALUE PROBLEM FOR PROPAGATION OF ELECTROMAGNETIC WAVES 

Calculations at 410 Me and 3000 Me using the bilinear Model have been continued, 
as described previously. 

Eigenvalues have been obtained for the Inverse-square Model at 50 Me using 
power series expansions. Work has been done toward slallar calculations using ssyMptotlc 
expansions and the calculation of normalized eigenfunctions by asymptotic and power series 
Methods. 

Programmers working on this problem are Professor H. B. Dwight and Dr. R. M. Ring. 

R. M. Ring 
Lincoln Laboratory 

199 N. STEADY LAMINAR FLOW OF A COMPRESSIBLE FLUID IN THE ENTRANCE REGION OF A TUBE 

Solutions for the first and second sets of the differential equations were 
obtained for the case where temperature dependence of the viscosity and thermal conduct¬ 
ivity of the compressible fluid Is taken Into consideration. 

More solutions are to be obtained for different entrance Mach numbers and 
thermal conditions at the tube wall. 

T. Y. Toong 

Mechanical Engineering 

219. TRANSPORTATION PROBLEM 

Changes have been made In the program to aolve the classical transportation 
problea by the stepping-stone method, (see Summary Report #45) which will allow the 
storage of coat data without specifically storing costs for shipping routes which are 


where G, H, and I arc known functions of the initial conditions and j is a unit impulse 
function. This program is nearly checked out. With this form of the equation the 
solution will be a Green's function which can be used In obtaining iterated solutions 
of the original equations. 

Preparations are being made to combine the various parts of the program which 
have been successfully checked out so that the solutions to the eauatlons can be computed 
in a continuous iterative process. 

Programmers working on this problem are P. Castillo, M. Jacobs and D. Cooley. 

D. R. Cooley 
Meteorology 

234 N. ATOMIC INTEGRAL 

A subprogram has been written and tested which evaluates the radial Integrals 
arising from matrix elements of 


d 1 



where r is distance measured from a point X Itself at distance R from the origin of 
coordinates. The radial integrals Involve products of two normalized atomic orbitals 



and the radial operator 
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The*. Integrals are evaluated by simple recurrence relations slnllar to those 
for the operator l/r , described elsewhere.* 

R. K. Resbet 

Solid State and Molecular Theory Group 


• Quarterly Progress Report, Solid State and Molecular Theory Group, October 15, 1955 
P. a 


236 C. TRANSIENT KESPON8E OF AIRCRAFT TO HEATING 

During the past quarter the WWI has been primarily used to evaluate character 
istlc values and characteristic vectors for matrix equations of the form Ax ■ \Bx. 

The program developed by F. J. Corbato and described In Memorandum DCL-58 has been 
employed. 


245 N. THEORY OF NEUTRON REACTIONS 


and 


taring this period the program was revised to give greater accuracy for ^ ■ 5 
d ■ 6, to make the cross sections a function of Instead of and to let 


5 ■ 


C 


where C is a constant. Several cases with x • 0 and C 
(Aj the basis of these a complete case was done for C ■ 
same C and 3 cases were done varying and letting 


4 . 

2 


42. 


and 5 as parameters were run. 
1.65 and J ■ ,08. For this 


From these,angular distributions were calculated in order to moke comparisons with 
experimental values. 


The work accomplished is reported In Wright Air Development Center Technical 
Report 56-287, entitled, Application of the Variational Method, the Galerkin Technique, 
and Normal Coordinates In a Transient Temperature Distribution Problem’, which is soon 
to be published. 

At this time, it appears that future use of the ill Computer on this problem 
will be limited to a) production runs of existing programs; b) determination of character¬ 
istic values and vectors for matrix equations; c) inversion and multiplication of matrices. 

L. A. Schmlt 

Aeronautical Engineering 

241 B. TRANSIENTS IN CONTINUOUS DISTILLATION SYSTEMS 


Future plans are to do more cases In the vicinity of the present C and _J 

E. Campbell 

Nuclear Science Laboratory 

246 B,N. SCATTERING FROM OXYGEN 

This problem is now terminated. The results can be found in a Ph.D. thesis 
entitled "Polarization and Exchange Effects in the Scattering of Electrons From Atonic 
Oxygen' submitted by Aaron Temkin to the MIT Physics Department, May, 1956. 

A. Temkin 
Physics Department 


A study of linear approximations in solution of control problems and problems 253 N. AN AUGMENTED PLANE WAVE METHOD FOR IRON 

involving the rate of change from one equilibrium state to another has been expanded. 

Correlations based on these approximations have been developed for determining the Little machine time has been used during the past quarter. That used has been 

approximate of a column as the time after a sudden change in operating variables becomes on a few calculations on the fine structure of atomic iron. 

large. 

Future plans are to obtain the energy band structure of face- and body-centered 

A final report on this problem should be completed shortly. iron. 


S. H. Davis, Jr. J. H. Wood 

Chemical Engineering Solid 8tate and Molecular Theory Oroup 

244 C. DATA REDUCTION FOR X-l FIRE CONTROL 256 C. WHIRLWIND I- UNIVAC SCIENTIFIC 1103 INPUT TRANSLATION PROGRAM 


This problem is concerned with computing from fire control signals how far 
fictitious projectiles miss an observed target. The program tape was ready prior to 
this quarter. Work during this quarter consisted In continuing to reduce new target 
run data by use of the program tape. 

Future plans for th* problem again call for .till acre data reduction, perhaps 
also using supplementary program, (y.t to be prepared) for clearer presentation of the 
reduced data. 

J. M. Stark 

Instrumentation Laboratory 


taring the past quarter, work described in Summary Report No. 45 on the WWI- 
1103 Input Translation Program to enlarge the vocabulary of the system has been completed. 
Numbers of the general form 

± 12.3456 • 10* • 2~ 3 

may be written where Integer number, only were previously allowed. 

A Digital Computer Laboratory report describing the entire system will be 
completed during the next quarter. 

J. M. Frankovlch 
Lincoln Laboratory 
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260 N. KNEHGY LEVELS OP DIATOMIC HYDRIDES 

Results of the work to date have been submitted as a thesis to the Physics 
Department at MIT. Por details see the Quarterly Progress Report, Solid State and 
Molecular Theory Oroup, MIT, July 15, 1956. p. 8. 

The problem Is being continued with the hope of calculating other molecular 
quantities, such as dipole moment. 

A. J. Preeman 

Solid State and Molecular Theory Group 

261 C. FOURIER SYNTHESIS FOR CRYSTAL STRUCTURES 

During the past quarter, Professor Buerger continued with the refinement of 
the structure of Wollastonlte. 

The structures of Diglyclnhydrochloride and -bromide are completely solved. 

A very accurate solution has been obtained for the chloride. Both structure determinations 
are now prepared for publication. 

Programmers working on this problem are Professor M. J. Buerger and T. Hahn. 

M. J. Buerger 
T. Hahn 
Geology 

262 N. EVALUATION OP TWO-CENTER INTEGRALS 

Some of the exchange and hybrid integrals which were evaluated by Merryman's 
program were found to be inaccurate. These Integrals have been recalculated by 
P. J. Corbato’s routine. The Whirlwind time per integral Averages about 80 seconds. 

H. A. Aghajanian 

Solid State and Molecular Theory Group 

264 C. OPTIMIZATION OP AIRCRAFT ALTERNATOR REGULATING SYSTEM 

The program for minimizing a function of n variables in the presence of equality 
constraints by means of the method of steepest descent has been rewritten in order to 
utilize the secondary storage units of the computer. Thus, it is now possible to solve 
problems of relatively large magnitude. Changes in the order of computations have been 
made to reduce the time required to achieve a solution. Thia new program has been tested 
and proven capable of handling a 27 variable problem. 

In the next quarter, the optimization of an aircraft alternator with respect 
to weight will be run. In addition, work on Improving the methods of computation will 
continue. 

Programmers working on this problem are R. R. Brown and J. B. Dennis. 

R. R. Brown 

Energy Conversion laboratory 
Electrical Engineering 
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266 A. CALCULATIONS FOR THE MIT REACTOR 

Por this past quarter, the most work has been done on improving the input 
data to refine answers. Future plans are to continue this improvement throughout the 

sumaer. 

T. Cantwell 
Chemical Engineering 

270 B. CRITICAL MASSES IN D O MODERATED REACTORS 

During the past quarter, the following programs were completed and production 
runs finished: 1) Serber Wilson sphere, 2) Fourier 3-group sphere, 3) Two-group 
cylinder, 4) Fourier 3-group (approximate) cylinder, 5) Fourier 3-group (exact) cylinder. 

The data is consistent with theory for the most part. 

Future plans include further debugging of the multi-group transport program 
for spherical reactors. 

J. K. Powell 
Chemical Engineering 

272 L. GENERAL RAYD18T SOLUTION 

This problem was solved satisfactorily by the method described in previous 
Summary Reports. A full description of the method and program can be found in Lincoln 
Laboratory Memorandum 2M-0508, "Solution of General Raydlat Equations on Whirlwind I". 

A. Zabludowsky 

Digital Computer Laboratory 

273 N. COSMIC RAY AIR SHOWERS 

In the last Summary Report we mentioned that we were going to change the 
empirical function that is used in the least-squares curve-fitting part of the program. 
This has been done and now there is a unique value of each parameter of the fit for a 
given fit. The resolution of the method we use will be determined by a Monte Carlo type 
of calculation in which artificial showers will be made up and Polsaon fluctuations will 
be put in. The prograa is virtually in its final form now. The computation tlae per 
shower is about 45 seconds and any changes in the prograa will probably be for the 
purpose of further reduction of the time. 

Most of the computer time used by us in the past quarter was devoted to 
analyzing air shower data. This will be the main effort in the next quarter. There is 
a backlog of over 1000 showers to be analyzed. 

P. Scherb 
Cosmic Ray Oroup 
Physics 

278 N. ENERGY LEVELS OP DIATOMIC MOLECULES (L1H) 

We are continuing with the self-consistent field molecular orbital treatment 
of the electronic energy. The binding energy for the ground state has been obtained for 
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fiv. Intwrmiclear dlatancea. Thw .ol.cular orbital approach, howavar, bwcowea lncraaa- 
lnaly poor aa the lntarnuclaar separation Increases, and It will be necesaary to apply 
additional configuration Interaction to the aolecular orbital ground state calculation. 
Work Is continuing at nore extended lntemuclear distances using the conventional con- 
figuration Interaction scheme dlacuaaed In previous reporta. 

Frogranaera working on thla problem are A. M. Karo and A. R. Olson. 

A. M. Karo 

Solid State and Molecular Theory Group 

288 N. ATOMIC 1AVK FUNCTIONS 

Calculations have continued on the hyperflne structure of lithium and boron, 
and on a wave function of high accuracy for the ground state of helium. 

Calculations on 4-electron molecules by the one-center method are still In 

progress. 

Approximate wave functions are being obtained for several configurations of 
Fe and for Mn, Mn**, and Cr*"*. 

R. K. Nesbet 
R. K. Vatson 

Solid State and Molecular Theory Group 

290 N. POLARIZATION EFFECTS IN THE FLUORINE ION 

The numerical integration of the equations for the perturbed functions V 
still show considerable sensitivity to the smoothness of the potential and the 
mesh size (we are using the Kutta-0111 library subroutine). However, for the uniform 
field case we have obtained a complete check on the perturbed functions of d-llke 
symmetry and considerable improvement in the p- and s-like ones. Recently we have 
obtained the V function of d-like symmetry for the dipole term of a point charge 
perturbing potential at a finite distance, R, from the origin. For R = 1.76 atomic units 
the V function for Ne has a behavior near the origin very similar to the corresponding 
uniform field (R «= ao ) case and a maximum at approximately r = R. 

Because of the delays and difficulties in obtaining our distorted functions we 
have undertaken a calculation for the HF molecule based on a small set of unperturbed 
functions to describe the F~ ion and including a d-llke distorted function similar to 
the one described above. To carry through this problem we have used R. K. Nesbet's 
tapes for the atomic integrals, potential Integrals, A matrix preparation, Roothaan 
scheme and integral transformation. These programs have all been described by R. K. 
Nesbet in previous Whirlwind I Summary Reports. 

L. C. Allen 

Solid 8tate and Molecular Theory Group 

293 C. ROLLING BEARINGS 

Progreaa h»» been held up while a new routine for non-linear equation* ia 
under inventlgation. If thla new routine la proved uaeful, it will be applied to the 
bearings equations. 
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297 C. DIFFUSION BOUNDARY LAYER 

The effects of mass transfer have been studied numerically in connection with 
40 analysis of the binary-mixture boundary layer. A system of equations describing mass, 
momentum, and energy conservation for the high-speed mixture boundary layer was numerically 
integrated to yield "similarity solutions for the case of relatively large and small 
densities for the injected mixture component. Beneficial alteration of the heat transfer 
r *te, skin friction, and ultimate equilibrium ( 'recovery") temperatures at supersonic 
velocities were shown to result upon employing low density coolants" such as, for 
example, helium. Primarily, the thermal effect is a consequence of the relatively large 
thermal capacity for such gases and a major requirement of the diffusion analysis is the 
determination of the specific heat variation through the layer. 

The Glll-Kutta procedure was used for integration purposes. Simultaneous 
solution of the differential equations (of a total order of 7) with two-point boundary 
conditions was completed with the aid of a linear influence coefficient routine. 

Successive integrations were programmed with initial estimates being self-corrected so 
as to yield proper end-point conditions. For a majority of the cases considered, 
approximation of the conditions was sufficiently reliable as to result in convergence 
within a time period on the order of a few computation periods of the results for the 
exact values. 

This problem is now terminated. A more complete description of the results will 
be found in a doctoral thesis submitted by the author to the MIT Aeronautical Engineering 
Department, June 1956. 

J. R. Baron 

Naval Supersoulc Laboratory 

300 L. TROPOSPHERIC PROPAGATION 

The last Summary Report on Problem 300 stated that tables of values of the 
Fresnel integrals for complex arguments had been computed using a convergent power series. 
These tables have been extended by using asymptotic series to regions in which convergence 
of the power series is too slow for machine computation. These tables are to be extended 
in the coming quarter. 

To simplify programming, closed subroutines using program parameters have been 
written and tested which carry out the common arithmetic operations on complex numbers. 

The use of these subroutines permits a program to be written in which real and immaglnary 
parts of complex quantities do not need to be considered separately. It is planned to 
construct interpretive subroutines for complex numbers, using a complex MRA so that a 
single operation nay be packed into one register. 

A closed subroutine has been written which obtains all the ordinates and 
abscissae of the extreme of a real single-valued function of a real variable throughout 
a specified range of the Independent variable. 

The autocorrelation function and power density spectrum of several samples 
taken from SHF signal strength records have been computed. Examination of the results 
indicates that Improved methods of statistical analysis must be employed. 

A number of matrices, ranging in order from six to twenty-five, were solved 
for their largest Eigenvalue and Eigenvector The results of these computations are to 
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A. J. Shaahaty 
Lubrication Laboratory 
Mechanical Engineering 
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b« used In a report on the design of scatter antennas. 

The LSR MA 1 was used to compute these Rlgenvaluts and vectors. An error 
In the coding of this LSR vas found, and has been corrected to give satisfactory results. 

Prograaaers working on this problea are W. Mason, P. IXiffy and J. Skenlan. 

V. Mason 

Lincoln Laboratory 

306 D. SPECTRAL ANALYSIS OP ATMOSPHERIC DATA 

The prograa for this problea has been successfully checked out and during the 
quarter, the routine use of this prograa was begun. 

It is anticipated that during the next quarter the extensive coaputations for 
which the prograa was written will be coapleted. 

B. Saltraan 

Meteorology 

309 B.N. PERFECT AND ALMOST PERFECT KC1 

For the perfect crystal the calculation of aatrix eleaents was coapleted, 
and the energy bands were calculated. The energy levels were calculated for twenty 
Independent values of the one-electron propagation constant k. Thus for each of these 
values of k an 8 x 8 secular equation with overlap terns vas solved, using a prograa 
written by F. J. Corbaco; this order results froa inclusion of eight different types 
of free ion orbitals, i.e., Cl“ 3px, 3py, 3px, 3ps, and K+ 3px, 3py, 3pr, and 3ps. The 
resultant Cl"3p band had a alnlaun at the Brillouin Zone Pace in the (111) direction in 
k space and a aaxlaua of very little curvature at the sane point. The aatrix eleaents 
involved were used to calculate a value for the cohesive energy of KC1; the result 
coapared well with experiments and with previous calculations done in the sane general 
Banner. 

For the Inperfect crystal the calculation of additional aatrix eleaents is 
almost coaplete and preliminary secular equations for energy levels have been solved. 

This work will be coapleted during the coning quarter. 

L. P. Howland 

Solid State and Molecular Theory Group 

312 L. ERROR ANALY8IS 

Ttoe following nodlficatIon has been made on this problea, as first described 
in Sunaary Report No. 43. A three dimensional syatea In space rather than a two 
dlaenslonal system on a surface leads to the solution of a 3 x 3 matrix (rather than 
a 2 x 2) for the eigenvalues and one eigenvector as final results. There will be 
other modifications and changes with the goal of studying different systems. 

Programmers working on this problea are L. Peterson and E. Hutcheson. 

L. Peterson 

Lincoln Laboratory 


WHIRLWIND CODING AND APPLICATIONS 


317 C. EXTRACTION OF STABILITY DERIVATIVE FROM FLIGHT TEST DATA 

The non-linear pitching moment and lift coefficients have been successfully 
extracted froa simulated flight records over an angle of attack range from -0.1 radian 
to 40.25 radians, A three second response and a time Increment of .1 second between 
data points were used. Programs are now being written which will permit the extraction 
of the derivatives over all angles of attack of Interest and also for the extraction of 
the lateral stability derivatives. 

A program written by Hiss Eatherlne Kavanagh of the MIT Dynamic Analysis and 
Control Laboratory has been used to obtain the simulated responses for the above studies 
and also has been used to verify analytic studies of the cross coupling between the 
lateral and longitudinal modes of a high-speed fighter aircraft. 

Programmers working on this problem are T. M. Carney, L. L. Maxxola, 

M. N. Springer, and L. E. Wilkie. 

L. L. Maxxola 

Aorophyslcs Research Group 

326 C. PRODUCTION FOR TRANSPORTATION PROBLEM 

Work done under problem 326 during the past quarter has been reported on In 
problea 219. 

W. S. Jewell 
Electrical Engineering 

327 L. PREDICTION ANALYSIS 

The problem , as described In Suastry Report No. 45, la In working order. 

There was one modification completed, with respect to alternate methods of computing 
certain functions. During the past quarter the major activity has been production runs. 
For the future, the main body of this problem will be Incorporated In a much larger 
systesM study. 

Programmers working on this problem are L. Peterson and K. Hutcheson. 

L. Peterson 
Lincoln Laboratory 

330 C. POSTFAILURE RESPONSE OF AIRCRAFT STRUCTURES 

Numerical solutions for the postfailure response of a beam represented by four 
lumped masses have been obtained using a finite difference solution of the equations of 
motion. Several simplifications have been explored and have been found to be In good 
agreement with the "exact" numerical solution. In addition, the WWI solution has shown 
satisfactory correlation with some experimental results. It Is contemplated that this 
study will be extended to Include the postfailure response of structures that are char¬ 
acterized by plate-1Ike deformation patterns. 

R. D'Aaato 

Aeroelastlc and Structures Research Lab. 
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333 A, BOUNDARY-UL YER CHARACTERISTIC3 OF A STEADY LAMINAR FLOW OF A COMBUSTIBLE MIXTURE 

OVER A MOT 8URFACE 

Solutions wsrs obtained for two specific values of the surface teaperature 
and for a second-order reaction of the Arrhenius type. 

The effect of activation energy is to be studied in the future. 

Prograaaers working on this problea are Professor T. Y. Toong and A. Shashaty. 

T. Y. Toong 

Mechanical Engineering 

334 C. PARAMETRIC STUDY OF MODAL COUPLING AND DAMPING 

Previously a program had been formulated which would yield the simultaneous 
solution of three coupled vibratory equations. This program gave the modal response 
in each of the three modes for one set of initial conditions. During this quarter, the 
program was altered to yield solutions for five sets of initial conditions in a single 
run. 
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improve the recurrence formula by using higher order approximations to the boundary 
conditions. 

S. H. Crandall 
Mechanical Engineering 

341 C. STATISTICAL AND DYNAMIC FORECASTING MXTHOU6 

As explained In the previous Summary Report, the major effort is being devoted 
toward the systematization of a method for multiple linear prediction. The programs 
for computing the covarlancea and performing the diagonallxation were completed. These 
programs were applied successfully to ses level pressures over a network of 64 stations. 
In the future it is expected thst these programs will be applied to other eeta of data, 
l.e. t temperature, pressure change, etc. 

This method was applied to the functions (Z's) developed by the Synoptic 
Climatology project. The functions chosen were those computed for the 700 mb 5-day 
mean height anomalies. This approach to 5-day prediction will be expanded by the 
incorporation of temperature parameters. A final expression for the prediction equation 
will then be formed. 


Two production runs, each for a specific airplane altitude, have been run 
successfully on a given airplane for five airplane velocities. Future plans are to run 
similar parametric studies on several other airplanes. 

K. Wetmore 

Aeroelastlc and Structures Research Lab. 

338 C. OPTIMIZATION OF RAM-AIR COOLING SYSTEM 


Processing of cllmmtologlcal data for multiple linear prediction by these 
techniques has demanded continued use of Whirlwind 1. 

Numerical experiments are being continued to find s satisfactory simplified 
meteorological forecasting model suitable for generating a sizeable climatic record 
through extended integration with respect to time. This climatic record may be par¬ 
ticularly useful in Investigating properties of least squares prediction as applied 
to a non-linear system similar to the atmosphere. 


Uirlng this quarter, work on problem 338 was completed. The results of the 
project, whlcj will be published in a Wright Air Development Center Technical Report 
entitled: "Cooling and Materials Investigations for Aircraft Generator Equipment", 
are summarized as follows: i) Optimum ram-alr cooling systems have been found for a 
specific set of flight conditions in the range between Mach 1 and Mach 4, and Sea Level 
and 100,000 feet; 2) Subroutines for finding the weights of sir-to-sir heat exchanges 
and sir compressors have been developed; 3) The variation of loss of an aircraft 
generator, as the machine weight is changed, was computed. 

R. Moroney 

Energy Conversion Group 

Electrical Engineering 


339 A. VIBRATING BEAM 

A recurrence formula for integrating 

d 2 1 

+ 5 1 - 0 


vs. exp.rlaent.lly found to be unstable. A program for finding the eigenvalue, of the 
recurrence formula wee written and run. The reaulta of thla explained the lnatablllty 
and Indicated how thw recurrence formula should be improved. It la hoped to further 


dx 4 


The work being done 1» under the .upervision of Profe.aor Edward N. Loren*, 
meteorology Department. 

Programmer, working on thl. problem arm B. A. Kelley, B. Shorr and K.Bryan. 

K. A. Kelley 
Meteorology 


342 B. TRANSIENT HEATING OF SOLIDS 

When a alab having a constant Initial teaperature Is suddenly brought Into 
contact with another medium having an Invariant teaperature (either more or less than 
that of the slab) a floe of heat will occur within the alab with accompanying temperature 
changes. The problea of what the teaperature Is at any apeclflc point demands knowledge 
of the period of heat flow, physical propwrtles of the material, the Initial slab 
temperature and surrounding medium's tempersture which are embodied in Fourier'e General 
Law of Heat Conduction, or 



While the boundary conditions specify 

dt 


k I =Xh(t - t_) 

* | x=*L s 


the resistance to heat flow, 
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This problem la solved using the simplified assumptions associated with semi- 
infinite and finite Isotropic, homogeneous solids with uni-dimensional heat flow, l.e., 
physical properties are constant Irrespective of temperature. 

Defined mathematically, a seml-lnflnlte solid Is bounded by planes at x » 0 
and z a as while a finite solid has boundaries at x a -L and x a +L. 


Without a surface resistance to heat flow (h a <x>, a = O), a solid Immediately 
upon contact with the surrounding medium assumes the external temperature upon Its faces 

so that t ■ t . 

s a 


airing the past quarter values were obtained for the following equations and 
also for the first root of m 3 a cotp. 


Y 

avg 



, -(21+1) X 

(-1) e _ 

(21 + 1 ) 

2 Tr 2 

-(21+1) 2 JL X 
( 21 + 1 ) 2 


Uirlng the next quarter It le hoped that certain errors In the Iteration 
routine for = cotp will be Ironed out so that the following equations can be readily 
solved In a Banner similar to that used for equations 1. and 2. 
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343 C. WEATHER PREDICTION 

Uirlng the current quarter the research prograa has been concerned with two 
phases of Weather Prediction. 

Data for December and January 1948-52, Inclusive, were utilized to classify 
weather maps Into three classes, according to the vertical structure of the weather 
systems. Linear prediction equations were developed for each class for the determination 
of pressure change. An interesting and unexpected result was obtained. It was found 
that the accuracy of the prediction differed markedly between classes. Attention Is now 
being directed toward an explanation of this result since It should lead to more accurate 
weather prediction by statistical methods. 

The research on the development of new predictors has been completed. With the 
new predictors it is now possible to work directly with data from reporting stations. 

It is no longer necessary to draw weather maps and read off grid-point values. Further, 
these new predictors make it possible to represent the maximum amount of variance of a 
field with a minimum number of predictors. 

In the future it Is planned to pursue the work on linear prediction for various 
classes of weather systems. Data is being processed at the National Weather Records 
Center in Asheville, N.C. for inclusion in this program. The results to date Indicate 
that this line of approach Is a fruitful one. 

Work on this problem has been under supervision of Professor James M. Austin, 
Meteorology Department. 

Programmers working on this problem are E. A. Kelley, B. Shorr. H. Brun and 
R. Huschke. 


x - 1 L Y 

8 


oo 

i-1 <m 2 f^ + ; 
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E. A. Kelley 
Meteorology 


345 B. MATRIX MULTIPLICATION 


x - 0 Y a 2 


1-1 (m 2 p 2 + ■ + Dcos^ 


< 2 




Y 3 2 2. -5=_ 

,vg 1-1 e* t" 2 0i + ■ + i] 


The portion of the problea which was used for a doctoral thesis Is now 
completed. The thesis was submitted by R. R. Archer to the MIT Mechanical Engineering 
Department. 

In the future, computations will be completed to prepare this problem for 
publication. 

Programmers working on this problem are D. Grlne and R. Archer. 

R. R. Archer 
Mechanical Engineering 


346 B. ANALYSIS OF COMPLEX SPECTRA 
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R. Gel man 

Mechanical Engineering 


The object of the analysis of the spectrum of an atom Is to determine the 
spacing and quantum numbers of the energy levels of the atom. In principle, this can be 
done by measuring the wave length of all the spectral lines emitted by the atom, convert¬ 
ing these save lengths to vacuum wave numbers, and then performing all possible subtractions 
of these wave numbers and looking for repeated differences, 

A complex spectrum such as Is emitted by the rare earth elements contains several 
thousand spectral lines. Since none of the wave lengths can be measured exactly, many 
repeated differences will be obtained when performing the subtractions which do not 
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correspond to true energy level separations In the atom. For this reason, additional 
data suet be obtained. These data are usually obtained by Beans of Zeeaan Spectra In 
which the sample is excited in a magnetic field. The field will cause the spectral 
lines to split into many components. The pattern thus obtained enables one to neasure 
two quantum numbers of the energy levels Involved in the transition. The consistency 
of those quantum numbers in the various transitions will separate the true energy level 
spaclngs from the false ones. 

The Whirlwind I computer mill be used in this problem to convert the measured 
wave lengths to vacuum wave numbers by solving a three term polynomial equation, perform 
all the subtractions, search for repeated differences, and check the conslst&ncy of the 
Zeeaan data. This will be done for Erbium which has approximately 2500 spectral lines 
in the visible and near-ultraviolet. Once the program has been set up, It can be used 
for the analysis of any other complex spectrum. 

IXjrlng the past quarter, the analysis of the spectra of two elements, Br and 
Pr, was started. The computer was used to search the Zeeman data of these spectra for 
repeated values. Then the difference in the wave numbers for these lines was computed. 

The energy differences which are needed to complete the analysis are found by searching 
for repeated wave number differences. None of these wore obtained in the initial analysis. 
This Is believed to be due to the fact that the data available at this time Is not com¬ 
plete enough to give the desired results. The computer was then used to search the 
complete list of wave numbers for these elements for selected repeated differences. 

Several repeated differences were found which appear to be significant. 

The next step in the analysis consists In obtaining additional Zeeman data 
and searching for any additional energy level differences. It now appears as if the 
doslred results will be obtained in the near future. 

J. W. Lindner 
Spectroscopy Laboratory 
Physics 

350 D. LONG-RANGE FORECASTING 

The 30-day forecasting problem is to be attached by means of a method for 
specifying (or reconstructing) monthly mean pressure anomaly maps in a completely object¬ 
ive mathematical way. The method Involves the use of empirical orthogonal functions”, 
determined from a long series of past data, which will perform analogously to (but more 
efficiently than) analytical orthogonal functions such as spherical harmonics. The 
coefficients of the empirical functions will be put Into least-squares linear regression 
forecasting equations. Aside from forecasting possibilities, this approach will give 
us much valuable information about the characteristics of long-period mean maps. 

The original description of the problem stated that wwx will compute as a 
first step a network of covariances; later it will be used to diagonalixe the matrix of 
these covariances--a matrix which may be as large as 100 x 100. This matrix has been 
converted to a 90 x 90 matrix of correlations. Because of its size, this matrix must 
be partitioned and diagonalized in sections, only the major eigenvectors and eigenvalues 
being retained. 

The three sets of variances and covariances required in the original statement 
of the problem nave been computed. These have been combined and converted by hand into 
a 90 x 90 matrix f correlations. The diagonallzation of partitioned sections of the 
matrix has beer compJ^ 'd. The recombination of the matrix sections and final dlagonal- 
lzatlon will occupy *he next few weeks. The set of eigenvectors so obtained will 


46 


WHIRLWIND CODING AND APPLICATIONS 


constitute the "empirical orthogonal functions" mentioned in the original problem state¬ 
ment. The determination by matrix multiplication of the amplitude coefficients of these 
functions is the next major task for machine computation. 

Programmers working on this problem are E. Kelley, B. Sborr, and D. L. Gilman. 

D. L. Gilman 

Meteorology 

351 B. NON-UNIFORM FUEL DISTRIBUTION 

A general computational program based on the finite differences approximation 
of simultaneous differential equations for a 2-group, multi-region reactor (Age-Diffusion 
equations) has been successfully set up. Fast and thermal distributions of modified 
fluxes ( r x Flux) for various critical spherical assemblies have been computed and 
checked with analytical solutions. 

In order to obtain s uniform volume power density in the core, the concentration 
of the fuel is deliberately non-uniformlsed in various zones of the core. The calculations 
are performed on s prototype high power reactor. IXie to very high flux level, fission 
poisoning will be taken into account. Criticality will be adjusted by varying the total 
critical mass and thickness or composition of reflectors. 

A. Butter 

Nuclear Engineering 

354 D. RESPONSE OF A SINGLE STORY CONCRETE BUILDING TO DYNAMIC LOADING 

The object is to study the dynamic response of the root of s single story 
reinforced concrete building for various loading functions. 

Tbe roof is a system of infinite degree of freedom consisting of a one way 
•lab supported on relatively stiff beams. The root system ess resolved Into • two degree 
ot freedom system by assuming deflection shapes for the one way slab and lt» supporting 
beams and using Lagrange's equation to obtain the two dynamic equations 


.. C 1 P + V. - Vi 

v * -r:- 


p - V. - Vb 


where y and y* are the accelerations of the center of tbe beam and slab raspectlvely 
* b • 

y b la tbe deflection of the center of the beam 

, 1. th. daflection of the center of the slab r.lativ. to the center of the b~m 

a 

P is ths total load on ths slsb area 
M 1 and M 2 are masses, and 
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V 


C , C , C , C are constants depending on the assumed deflected shapes and 
2 3 4 5 

resistance characteristics of the systea. 


357 B. PROPAGATION OF ROUNDOFF ERROR 
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The nuaerical integration portion of the prograa calculates the deflections 
of the beaa and slab for successive time Intervals by utilising the backward difference 
equation 


y (tfD * 2 y u> - y (t-n +y( *° 2 

Proa tbeae delfectione new resistances are computed and the process repeated. 
This nuaerical procedure Is repeated approximately 100 times per trial loading. 

The loading portion o l the program Is an empirical curve which la a function 
of the shock wave characteristics (i.e., square save, triangular save, etc. and shock 
strength). 


The control portion of the program consists of choosing a new shock strength 
for the next trial loading on the basis of the maximum deflections of the slab previously 
obtained. 


IXirlng the past quarter, the program has been entirely written and debugged 
and Is now running. A great many more runs will be made this summer with smaller pro¬ 
grams which are segments of the program which is now running. 

B. Landry 
Civil Engineering 


3SS B. QUANTIZATION ERROR 

This problem was proposed by Mr. Bernard Wldrow In his thesis, "A study of 
Rough Amplitude ftiantizatlon by Means of Nyquist Rampling Theory.'' He developed there¬ 
in a theoretical proof that quantisation error Is statistically random and Independent 
of quantizer input when Nyquest'a restriction, applied to amplitude sampling (quantization) 
Is followed, i.e., sampling rate should be at least twice as high as the highest fre¬ 
quency caaiponent of the probability distribution. 

Experimental check on this theory can be obtained by determining the effect of 
roundoff on the autocorrelation function of a set of data. This can easily be accomplished 
by the wwi computer, using the autocorrelation prograa devised by Mr. Douglas T. Ross 
of the MIT Servomechanisms Laboratory, and a supplementary program causing repeated 
determination of this function (about five times) after binary roundoff, one bit at a time. 
The critical sampling rate thus determined can then be compared with that predicted by 
the Nyquist criterion, as a measure of the theory above. 

This problem was concluded after determination of the autocorrelation functions 
of error due to roundoff for three sets of data. These error functions, deduced from 
computed functions of the data before and after roundoff, agreed with Mr. Wldrow's theory 
that this error Is uncorrelated and of flattop probability distribution. The results 
of this problem were submitted a. a Bachelor's thesis to the MIT Department of Electrical 


The flow graph for the second order difference equation to be studied la 
shown in Fig. 1. At each multiplication (i.e., ...6, 1.2, -.8), artificial quantisation 
la introduced by Beans of the shift orders. The Input supplies random numbers, in order 
to obtain statistics for the output. An unquantlzed solution is obtained for each step, 
in addition to the quantized solution. The difference between the unauantlzed and the 
quantized solution Is foraed. The mean and second-moment for this difference is obtained 
for many different runs at each of the 50 steps of the solution. Also, a histogram of 
the 50th step error la kept. 

The problem will be solved, using different levels of quantization. The 
experimental results will be compared with the theoretical predictions of B. Wldrow's 
Doctoral Thesis Proposal, A Study of Rough Amplitude Quantization by Means of Nyquist 
Sampling Theory ' (Memorandum 6M-3862, Lincoln Laboratory, MIT), 



Fig. 1 Second Order Difference Equation 


The problem was completed during this quarter. The results compared favor¬ 
ably with theoretical predictions based on B. Wldrow's doctoral theala. The results of 
the problem, together with the analysis, appeared in my bachelor's thesis (Propagation 
of Quantization Error ’, S.B. Thesis, MIT, Department of Electrical Engineering, June, 
1956). Since the desired results were obtained, there are no future plana for the 
problem. 


A. I. Green 
Electrical Engineering 

358 B. VERTICAL TAIL LOADS DUE TO ROLLING PULL-UP 

The purpose ol ttale problem le to find the loads In the vertical tall of my 
16 72 airplane design. The critical maneuver to be Investigated aay be described as a 
rapid roll in the presence of a high normal load factor. The maximum tall load obtained 
for thla rolling pull-up will then be compared with tbet value found by neglecting non¬ 
linear and quasi-linear terms of the system. Bron this comparison we should be able to 
determine how large an effect the non-llnearlties of the system have on the vertical tall 
load. 


H. Pople 

Electrical Engineering 


The Lateral equations of motion of the airplane are: 
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.00263 + .0322 <P + .00101 -.02SS*f> + .0260 - .002296^(0 

.00101 if - .00991 (P + .0115 if + .01167/7 - .0630 » -.OOOGSB^ft) 

-.272 elnf* + 2.037/> + .40 + *.030 - 0 

Applying the trapezoidal rule we have In finite difference form: 

£.00263 + .0322 <p „ + ^00101 - .0253 + ,026 j" " .°° 229 b^lt) 

-“"tPn-l^fft-lV •““&-! + 2^n-l]- O28 f 0 n-l + ^n-l + T K-l] 

£00101 - .00991 + [ 011S + - 116 iJV n * - 063 r ^n “ '• 00063 »» (t) 

+ 00991 [<Pn-l + ! ^n-l] -• 0116 [^ > n-l+ I ^n J+ 063 [ Vl H *»-l + T Vl] 

2 -° 3 2 + [•« |* + 3 03 aK * • 272 • ln f n -r a 03 & , o -l + i^n-J 

-•< fco-1 + *.-l + - 2 * 03 [’.-I + 2 *»-l] 

Then the tall load: 

L ¥ » a<f - b 0 

The aileron angle will lncreaae linearly from zero to maximum detraction In two-tenths 
of a second; from then on It will remain constant. The program will be stopped after 
the peak tall load is reached. Por the urst run an increment of .01 second will be 
used at early timo and a coarser Increment of .1 second at Inter time. There will be 
approximately 200 cycles. 

The work on this problem has been successfully completed during the past 
quarter. A complete description of this problem has been presented In a S.B. thesis 
submitted to the Department of Aeronautical Engineering. 

Nomenclature 

T • angle of bank 

"V = angle of yaw 

0 • sideslip angle 

8 a • aileron angle 

t = time 

h - increment of time 

L^ • vertical tall load 
a, b = constants 
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359 B. COMPUTATION OF STRESSES IN A TRANSVERSE WEB PRAME 

A transverse web frame with non-symmetrical loading can be approximated by a 
closed ring of variable cross-section. A closed ring Is a structure of three degrees 
of lndetermlnancies. A web frame has In addition two vertical pin ended stanchion 
supports which complicate the problem by introducing 2 additional degrees of lndeter- 
mlnancy . 

The ring was divided Into 49 divisions of varying length. Loads (water pressure 
loads, weight loads and shear loads! were converted Into horizontal and vertical components 
at each station. Moments at each station were computed and Bruhn's adaptation of 
Castlgllano's method of least work was used to solve for the 5 redundante. Stresses due 
to bending and direct loads were computed and the process cycled to produce allowable 
values of these stresses at each station. 

• to - x ,) Zl y = (y - y ) + 25Ty^ 

n n-i n n -1 


sin 9 tdj/ds 


cos 9 =Ai/ds 


(sin 9) + (sin 9) 

n-i _n 


P = y X h x A n (Force due eater pressure) 


V « P (cos 9) H « P (sin 9) 

w w av w w av 


(cos 9! , + (cos 9) 

n-1 a 


M(about shell) 

N,A “ A(area of section) 


I„ . = I (about shell) -Ad 
N.A N,A 


M , P , Q . R and S are values of redundants. 
o o o 

u .p+p-Q+R +S+'S.'A 

(at each station) o oy ox xt x 1 


“l * “. + "L - “a - "v 1 “.(.hear) * 

* £<">„-! d Vl + 2 - (V> n-l d Vl 


Redundant unknowns were solved using Bruhn's eqatlons for conditions of slope 
and deflection. 

Bruhn's Equations 


triiuui d “H “ 

J * ds « 0 {Relative change In slope « o] 

o 

J'f y ds = 0 Reflection In x-dlrectlon - oj 
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H. Parechanlan 
Aeronautical Engineering 
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« y ds • 0 ritoflectlc 

o 

r*2 

J (i-*,)|d,.o [v 

'C (x ' V T d * = 0 [ Ve 


Deflection In y-dlrection 


Vertical deflection between B 


1«2 ’ °] 


Vertical deflection between B_ +• B 
• 3 




Method of Integratio n 


Trapezoidal Rule = 


< T -i da 

2_ [(Ord.) n + (Ord.) n J -5 


The above Integrations formed a matrix of 6 columns and 5 rows. The matrix 
was solved by Croat's method. After solution the redundants were Incorporated Into the 
existing horizontal and vertical loads. 

Direct loads » P « t-H cos 9 + £.V sin 9 

Shear loads = Q = S.H sin 9 + 2.V cos 9 

where £ H * P + 1 H i.V<s-Q + R + S +< V + w 

o w o w t 

Direct Btress « P/A Bending stress - My/I 
Shear stress = Q/A 
<T b = P/A + My/I 

A rational engineering stress allowance was used to compare with computed 
stresses and the section altered to conform to allowable stress. 

Miring this past quarter, the problen was completed and satisfactory results 

We ^.° bt ! i c ed- na r ® 8UltS °* the proble " * in b « found in a Naval Engineering thesis 
submitted by Brandt, Snyder and Thompson to the Department of Naval Architecture. The 
title of the thesis is Applications of Digital Computers to Naval Architecture Problem 

C. R. Brandt 
J. C. Snyder 
C. R. Thompson 
Naval Architecture 

303 A. EFFECT OF PRESSURE ON THE BENDING OF CURVED TUBES 

The bending of a thin walled toroidal shell under Internal pressure Is 
governed by the equations pressure is 

T" - qT + xR - 0 
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and the boundary conditions 

T(O) - 0 

S' (0) » 0 

Numerical solutions have been obtained by Clark and Relssner (Advances in Applied 
Mechanics II, Academic Press, N.Y., 1951, pp.93-122.) for the case q - 0 which 
corresponds to zero internal pressure, W C plan to obtain numerical solutions for 
several values of the pressure parameter, q. An iterative method based on a finite 
difference approximation will be employed, starting from Initial trial values of 
T (0) and S(0) the solution will be marched out and its asymptotic behavior compared 
with the desired behavior. Errors hero will be used to determine new Initial trial 
values for T (0) and S(0). Preliminary estimates indicate that adequate accuracy will 
be provided by an 0(h z ) process using 100 subdivisions or less. 

A program has been written and corrected and has been used to calculate T and 
2 for the case q a 0, The results of this calculation check other numerical solutions 
(for q = 0) available In the literature. 

Plans for the future are to use the established program to calculate T and 
S for cases where q / 0. 

Programmers working on this problem are Professor S. H. Crandall and Professor 
N. C. Dahl. 

N. C. Dahl 

Mechanical Engineering 

364 C. BLAST RESPONSE OF ROTOR BLADES 

This problem involves the determination of the stresses and the flapping of 
the blades of a helicopter rotor which while in flight is enveloped by a blast wave. 

The c.g. of the helicopter Is permitted to move In a direction parallel to the rotor 
shaft; non-uniform mass and stiffness distributions of the flexible rotor blades are 
taken Into account. In the Interest of simplicity, quasi-steady rather than unsteady 
airforces due to the blast are used In the analysis. Later the Influence on unsteady 
airforces on the response of the rotor will be studied. 

The transient rotor response is obtained by solving, by numerical Integration, 
nm + n + 1 simultaneous differential equations having variable coefficients, where m Is 
the number of blades on the rotor and n Is the number of normal modes of each rotating 
blade. The unknowns are the nm nodal amplitudes, n p-flapping motions, and the vertical 
translation of the c.g. By taking advantage of previous studies on blast response of 
fixed-wing aircraft, a somewhat simpler formulation of the above problen can be 
accomplished. This Is done by neglecting mass coupling and airforce coupling terms between 
the rigid-body response and the flapping and/or the nodal vibrations, as previous 
studies have shown to be permissible. In this way the vertical translation equation 
becomes a slngle-degree-of-freedom equation with the blast forces acting as the forcing 
function - it does not involve the flapping equations or the modal-vibrations equations. 

The flapping equations involve the flapping response and tho previously-computed rigid- 
body response; they are then solved. To compute the modal vibration, ths previously- 
computed rigid-body and flapping responses are Integrated. 

Whirlwind I is used to solve for the normal modes, shapes and frequmicles of 
the hinged, rotating, non-uniform (mass, stiffness, chord) blade. Whirlwind *s also 
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used to numerically Integrate the equations of notion; Integrations with respect to 
blade span and with respect to aximuttawise position of the rotating blade are performed. 

During the past quarter, programs have been written to compute the difference 
equations and matricos required to find the eigenvalues and eigenvectors which will be 
used to calculate the normal nodes and frequencies desired. 


The first phase of this problea will be completed sometime in the next 


quarter. 


This problea was initiated by J. Dugundl and K. A. Poss of the Aeronautical 
Knglneerlng Department, HIT. 

M. E, Callaghan 

Digital Computer Laboratory 


two solutions which are then converged, first, a cylindrical reactor with only a radial 
reflector la assumed and the exact solutions are solved. Second, s cylinder with only 
end reflectors is assumed, and a second set of exact solutions are obtained. The two 
sets of solutions are then made to converge. 

This problem dlfferes from number 270-43 because of the low enrichment in the 
uranium used as fuel. Allowance is made for the presonce of u23S by putting Into the 
program fictitious values of the constants for U^* 5 which are really effective values for 
the U235 - i;238 mixture. 

The following equations will be used: 

Tbese four boundary conditions result In 4 equations in four unknowns. 

By Cramer's rule, for non-trivial solutions, the determinant of the coefficients 
must equal zero. 


365. FINDING EIGENVALUES OF AN ASYMMETRIC MATRIX 

This problea Is concerned with finding eigenvalues of an asymmetric matrix. 

We will be using LSR MA 1 which finds the largest and smallest eigenvalues. A routine 
will be written to find the characteristic equation of a matrix using Frame's method. 

The characteristic equation will be solved by Hitchcock's method if possible. In 
particular, we have four 9x8 matrices to solve. 

Ikirlng the past quarter. Frame's method has been written and debugged, and 
does find the characteristic equations of the 9x9 matrices. However, we are having 
trouble solving the characteristic equations using Hitchcock's method. 

Ikirlng the next quarter, we will try to find a better method of finding eigen¬ 
values of asymmetric matrices In factoring polynomials. 

This problem is being done in close connection with the MIT Office of Statistical 
Services, whose Director, Dr. F. M. Verzuh, Initiated this problem on Whirlwind I. 

Richard Steinberg and Eleanor Donovan are doing the programming on the OSS IBM 650. 

J. Roseman 

Digital Computer Laboratory 

367 B. DETERMINATION OF CRITICAL MASS IN A CYLINDRICAL REACTOR 

The problem Is to apply two-group diffusion theory to the prediction of 
critical masses for finite cylindrical reactors which are completely reflected. 


The procedure is to write the diffusion equation 

^ *lj + “t) 0 - 0 


for the group of thermal neutrons and the group of fast neutrons. The resulting two 
equations containing average values of the constants tor each group are solved and 
appropriate boundary conditions applied. 

For the completely reflected cylinder being considered, exact solutions are 
not possible because of the inconsistency of the boundary conditions at the junction of 
side and top or bottom reflectors. Instead, a method of superposition Is used to give 


This determinant will be solved by Whirlwind I. 
End Reflectors Only 
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W - Cos (p^a) 


U - Slab (V z) 


V = Cosh (V s) 

z 


Radial Reflector Onl) 


Coro Flux: 


p 2 - slnh J2* [r + \ - •] 

S ■ 8inh ) 2, j/ +T 'x -■] 

S * 81nh ^2x [j + T . -*] 


Q 2 - Slnh 


Reflector Flux: 


* lc • AX + CyJ Cos (u^x) 
Slos: « 2c - S^AX + S 2 AyJ Cos (p b *) 


D | 0 l C - . D | “i 

lc ^r lr d r 


°2r "5T 


WHIRLWIND CODING AND APPLICATIONS 


Neutron Current Continuous 


Applying the sane boundary conditions as before but in the z direction; 

at z «■ ± H/2: 0 a 0 and 0„ . 0 

lc lr 2c 2r 


Fast: 0 » Fz Cos (p s) 

lr 1 » 

Sio.: 0 2r - js 3 FXj + Ox 2 J Cos (p B *) 

X a j (p r) 

o r 

Y - I <v r) r 

I 3 i (R + T >1 

*1 •». <-> - r i.Iri j 

» i r ^ o, <»+t)7 

Z 2 ' *o ( ^2r r) ' K o ? J 2r (R + T)j K o ( ^2r r) 


2 2 


- V 2 = V 3 +V * 


1 2 r 2 . 2 

' lr ‘ l r + •*, 


r 2 2 
K 2r + I*. 


Boundary Conditions: 


A * r ‘ V *lc “ 0 lr and *2c ’ ®2r ,lu * Continual. 


plus the additional condition that the flux equal zero at the outer surface of all 
reflectors, results in eight equations in eight unknowns. Again Cramer’s rule is 
applied to secure a non-trivial solution. 

For superposition to be valid both determinants must be satisfied simultaneously 
for the reactor to be critical. 

The solution procedure is to fix all the properties of the materials and the 
reactor dimensions. A critical mass is assumed, and the values of total buckling (p2 
and V 2 ) are computed. Iteration continues until the 4th order determinant is solved 
giving values of and •%) J. With these values and can be obtained. Finally, 
all these values are suostltuted Into the 8th order determinant. If it does not equal 
zero, a new critical mass is assumed and the calculations repeated until a predetermined 
limit of convergence is reached. 

The following symbols not elsewhere identified pre: 

_ _ .th . th 

0^ * neutron flux i energy group, J region 

« diffusion coefficient 

r * distance in radial direction 


diotance in longitudinal direction 
longitudinal buckling 
radial buckling 


Imaginary buckllr.gs 


S . S , S ■ constants which are coupling coefficients between fast and slow 
1 2 3 flu. 
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1 1 
*2 “*«• 

I _1_ 

2 = 2 

L 2r (thermal diffusion length) 


During the past quarter, a series of runs using three different programs have 
been made using data from an operational reactor to determine the success of each type 
of calculation in predicting the critical mass. As a rer'lt of these comparisons one 
program has been discarded. One of the remaining two Is undergoing corrections and 
modifications. The third program tried has given very accurate results to compare with 
the measured critical mass. 

Plans are to adopt the most satisfactory of the two remaining programs and to 
determine through Its use the critical masses of several types of reactor cores over a 
range of uranium enrichment up to a maximum of 20% U 335 . These results should permit 
selection of the optimum design. 

Programmers working on this problem are J. W. .. .mett, M. Troost and 
J. R. Powell. 

J. W. Barnett 
Nuclear Engineering 

369. TEMPERATURE DISTRIBUTION IN A BEAM 

Tills problem is concerned with the solution of a set of 30 non-linear equations 
which arose In connection with a study of boat transfer in a cylindrical ring. First 
order finite difference approximations were used to solve partial differential equations 
which were present in this problem. 

As part of this study the same problem was solved on the IBM Type 650 Magnetic 
Drum Calculator where the drum calculator was used essentially as a desk calculator to 
obtain check solutions for the Whirlwind program. That is to say, only a limited portion 
of the overall problem was solved — time levels 0, 1, 40 and 41 respectively. Time 
levels 40 and 41 were selected to overcome a common pitfall encountered preparing test 
solutions on a large-scale calculator, namely, that of using as a check-point values 
which are non-zero. Specifically, usual techniques for initiation of check solutions 
at the beginning of a problem are dangerous because a number of the values are zero or 
non-zero. Obviously, a zero result obtained by any correct or incorrect method gives 
no measure of the correctness of the program. Selecting step 41 eliminated the non-zero 
limitation. 


Results obtained on Whirlwind I and the 650 agree to 6 decimal digits which Is 
quite good considering the differences in the two interpretive systems used on Whirlwind 
and tho 650 respectively. 

Programmers working on this problem are A. Z&bludowsky of the Digital Computer 
Laboratory, and M. Hermann of the Office of Statistical Services. 

F. M. Verzuh 

Office of Statistical Services 


371 L. ATMOSPHERIC PROPAGATION OF RADIO WAVES 

The problem consists of the Investigation of the influence »f the dielectric 
inhomogeneities of the lower atmosphere and ionosphere upon the propagation of electro¬ 
magnetic waves to great distances and elevations. In the first phase of this work the 
stability of rays at low angles to the horlzou will be investigated from solutions of 
the non-linear differential equation: 


dh 

dR 


1 

n(h) 



r n cos 6 I 

o o _o J 

n(h) <r+h)'J 


2 



- 0 


(1) 


where R t ct, c is the velocity of radiation in a vucuum, t Is the time required for a 
disturbance to reach a point (R,h) on a ray trajectory whose angle with the horizon at 
(0,0) is h the height of a point (K,h) on the trajectory and n(h) is the refract¬ 

ive index of the atmosphere given us a function of height. r y Is the radius of the earth. 
Equation (1) is to be solved for a large number of functions n(h) corresponding to measure¬ 
ments of refractive index in the lower atmosphere and ionosphere at various wavelengths. 
These solutions are expected to provide important information about radio optics in tho 
atmosphere. Seasonal daily and hourly fluctuations in the trajectories about mean tra¬ 
jectories will be obtained. 


In the second phase of this problem, the astronomical refraction R of radio 
•tars at various wavelengths will be computed from the Integral 


R<€ ) 

o 




dh 


Fluctuations of R<€ ) will be obtained from the use of measured functions n(h). 
Additional phases of°this work are to be specified. 

Two methods for solving the refraction equation have been tested, — tho Kutta- 
Oill method and the Tchebyshev method of Integration. The latter method was found to be 
more rapid than the former and appears to provide a satisfactory degree of accuracy. The 
radio range versus height for elevation angles from 0° to 5° have been determined by this 
method up to heights of 300,000 feet for different seasons of the year. 


In the future we plan to: 1) extend the problem using the above method from 
heights of 300,000 feet on up througn the Ionosphere at a number of frequencies; 2) 
investigate, using the above method, the hourly variations in atmospheric refraction. 


W. Mason 

Lincoln Laboratory 
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372 B. DESIGN OF SPHERICAL SHELL SEGMENTS 

The problem proposed Tor solution on Whirlwind I Is part of a doctoral thesis, 
carried out under the supervision of Professor C. H. Norris, on the subject: The design 

of precast doses". The general Idea is to construct a reinforced concrete spherical 
shell out of Individual pie-shaped segments, which are assembled by means of prestressing 
to form the complete dome. Each such precast segment will be hoisted Into Its final 
position and be temporarily supported at Its crown and bottom ring. 

The present problem, to be solved on Whirlwind, deals with the determination 
of the membrane stresses In this precast spherical trianglo during construction. 


The range of 0, the angle from the vertical axis of rotation should be o<0<^ and 
values should be determined In intervals close enough to give an accurate plot of the 
function. 

It Is obvious , that In view of the very high values of n numerical integration 
must be used. 

However, this expression for the neridional stress given above still is a 
function of two constants of integration. To determine them the equilibrium conditions 
for the complete shell segment are used. In writing these out, one comes across 
expressions of the following form: 


In shell thoery it is normally found, that an explicit solution in closed form 
for the stresses in tho shell can be obtained only for the simplest forms of shape and 
loading. However, for the given shape of a spherical triangle, where there is no rotat¬ 
ional symmetry, one must use series substitutions and numerical methods to evaluate the 
stresses. 


The proposed method is based on the following theory: our object is to 
determine the membrane stresses in a shell segment under its own load, when supported 
at its crown and bottom ring. This state of stresses will be the same as in a complete 
dome in which the vertical planes bordering the individual segment form planes of anti¬ 
symmetry. For such a case of antisymmetric load, one can express the load and stresses 
in terms of a Pourler series and can thus integrate the equilibrium equations. As a 
result of these expressions the shear stress, for example, can be expressed in the 
following form: 
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This consists for each value of n of three terms, 
constant of integration and one a numerical value, 
evaluated for the following values of n : 


two of them proportional to a 
Expressions of this form must be 


Number of Segments 

k - 2 
k = 4 
k * 12 

k - 24 
k - 36 


18, 
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42 
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sin 0 d0 


N 00n CO “ 0 d0 


"m* co ‘ 0 <l 


Cos 0 ) d 0 


which again must be evaluated by numerical procedures for each value of n as detailed 
above. 


As a result of the integrations outlined ubove a series of curves will be 
plotted which will finally give all the necessary data for the determination of the 
membrane stresses in any segment of a spherical shell for any opening angle 0. 

During the past quarter, part of the program was developed and checked lor a 
particular case. Hie result obtained verified the numerical solution. 

In the future, the program will be set up for the determination of stresses 
in a more general form and with the results from the computer,graphs will be plotted to 
determine the membrane stresses at any segment of a spherical shell. 

Programmers working on the problem are E. Traum and Assistant Professor S. Nanyet, 

S. Namyet 

Civil Engineering 

373 B. FLUX LEVELING IN HOMOGENEOUS REACTOR - PART I 

Spatial distribution of neutrons in a nuclear reactor usually prevents full 
exploitation of the outer regions of the reactor, since heat generation, which follows 
the spatial distribution of neutrons, is limited by heat transfer and metallurgical 
problems at tho reactor center. Neutron distribution may be made more nearly uniform 
by decreasing moderation in the core and allowing moderation in the reflector to increase 
neutron concentration In the outer core. 

The diffusion equation, which describes neutron distribution, is solved for a 
spherical (one-dimensional) reflected reactor for the time independent case. Coefficients 
of this equation, which vary with neutron energy, are replaced with averaged constant 
values by considering n energy groups of neutrons, each of which thon has its own 
diffusion equation. 
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These n diffusion equations are non-dimensionallzed, converted to a standard 

* 8 ‘“ d “ rd lteratlon “««>d> “d integrated. WWI solves n equation, of 
the following fora by trapezoidal integration: 


Y(x) « exp( 




where x Is dlnensionless space paraneter r/R 
a 1s a aaterlal property of the core 

Z Is the source tens, representing fission neutrons, or neutron, aoderated 
to a lover energy group 

f Lvra lnciude - 

hand computations'(not T ’'° nt5 ’ h0UrS ,re re 0 ulr «<> *«r 

seconds, both " ,I 19 

Program. S'S ft?"™T;., l !TZct <13 rr e " J t. debug the 

previously been hand-calculated! aolutlon to the problem which had 

At present data are being collected to feed WWI for final results. 

R. V. Kennedy 
Nuclear Engineering 

377 L. COVERAGE ANALYSIS 

Of the ob.ervatIona 1 for^.eumed'radarConfigurations^ T*' “ d «« 

first part of the problem trigon^tr^-i^K J “ d tllctlcal wltuatlons. In the 
atlon of parameters. The program 1. desirafd to JuJ^ 10 " 8 *” 801 V9d toT each oombln- 
Each combination results in six associated! ! ? * ">«imum of 4096 combinations, 

drum for accuracy requirements in order to IT ! ValueB thEt Wl11 be stored on the 

r 0 i “th b :.r cod fro ° doubie — « to ----n“tr!d" th v :“ a “!a!! or ttg!.t.r e 

functional T values°to T* 1 " “ 8l - 

observations. ^ and characteristics of the distribution of the 

been compldd Pr ° ble " 18 ^ ^ Pr ° CeM ° f *•*»« Programmed. Part of the first half ha. 

Programmer, .orbing on this problem are K. Robeson and L. Peterson. 

1*. Peterson 
Lincoln Laboratory 


APPENDIX 


1. SYSTEMS ENGINEERING 

lnterventlod! 1 d', the “T Servo « ch ‘ nl “» Laboratory ha. been using the WWI manual 
ta o!er t y 6 ^ 1P “ nt ln the dev «l°P»ont of high speed data reduction 

0 d r tOT the ” to e *P“ d l»eir research Into computer applications, it 

tedrd dt/ e T tUe “ nU81 lnPUtS be nade on the WWI computer. 

r itl0n81 0n '°“ 8wltche8 ’ ~n» of the new programs will be so com- 
tlU to !! d t*,t° T’ 3 ' par “ eters th8t *»e only reliable way to Instruct thorn 

to have thid.™ d * d88lgned ■“•“onlc languages and translation programs. In order 
!!! ♦ eoneral language structure available on a manual Intervention basts, it Is 

necessary to have a keyboard such as a Flexowriter for direct Input to the computer. 

The LIT Scientific and Engineering Computations Group have contemplated the 
following applications for the new facility: 

1. Demonstration programs would be a great deal more effective if this form 
of input were available for control purposes. 

2. Typewriter Input for Comprehensive System Flexowriter and post-mortem re¬ 
quest tapes. Short program modifications and post-mortem requests can presently be In¬ 
serted ln the Insertion registers. However, errors are easily made because the re¬ 
quired vocabulary Is awkward. A typewriter input facility would make available a 
normal mnemonic vocabulary for such purposes. 

3. Experimental use of a typewriter facility for direct operator control of 
the computer. Here we would consider using the typewriter to replace the button- 
pushings required of the operator during normal operations. Vocabulary similar to that 
of director tapes and performance requests would be devised for these purposes. This 
could easily prove to be an extremely convenient and efficient method of computer opera¬ 
tion. 

The new Input installation will be available for use by 4 July 1956. Much of 
the information to be Inserted via the keyboard will be the same as Is now Introduced 
via a free running photoelectric tape reader using punched paper tapes. The keyboard 
Input will also be treated as a free running device, i.e., selection of the facility 
by the computer may be followed by an arbitrary number of read Instructions, each of 
whlcb reads the next character which has been struck on the keyboard. The total equip¬ 
ment requirements amount to 15 relays and 20 tubes. 

WWI RELIABILITY 

The following is the WWI Computer reliability for the past quarter; 


Total Computer Operating Time ln Hours 
Total Time Lost ln Hours 
Percentage Operating Time Usable 

Average Uninterrupted Operating Time 
Between Failure Incidents ln Hours 

Total Number of Failure Incidents 

Failure Incidents per 24-hour Day 

Average Lost Time Per Incident ln Minutes 

Average Preventive Maintenance Time Per Day In Hours 
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2. VISITORS 


3. ACADEMIC 


Introduction 


Tours of the Whirlwind I Installation Include a shoeing of the flla, Making 
Electrons Count", a coaputer demonstration, and an Informal discussion of the aajor 
coaputer components. During the past quarter, 10 groups totalling 190 people visited 
the coaputer installation, representing the following: 


There are a number of graduate subjects In autosiatlc computation, numerical 
analysis, and now electronic data processing, offered at M.I.T. The present list of 
subjects directly related to machine computation Includes the following: 


April 23 
April 27 
May 3 
May 3 

May 8 
May 16 
May 17 
May 17 
May 29 
May 31 
June 8 


University of New Hampshire, Student Branch of the AIEE 
United Nations 

MIT Course entitled "Electronic Coaputational Laboratory" 

Class from Northeastern University 

U. S. Representatives of S.A.C.M , France 

Foxboro Corp 

MIT Class In Aeronautical Engineering 
MIT Class In "Machine-Aided Analysis" 

Gillette Co. 

Phillips Eindhoven, The Netherlands 
British Joint Services Mission 


The procedure of holding Open House at the Laboratory on the first and third 
Tuesday of each month has continued. A total of 155 people attended the four Open 
House tours during the quarter, representing aeabers and friends of the MIT Community, 
Weston High School, Hood Rubber, Northeastern University, Purdue, Tufts, the Massachusetts 
Department of Commerce, Dunlop of France, Godfrey L. Cabot Co., Rand, Inc., Mlcroflne 
Instrument Co., Massachusetts General Hospital, Edvard Devotion School, Raytheon Mfg. 

Co., Westlnghouse Electric Co., Remington Rand, Production Systeas, Inc., Elliot Bros. 
London office, General Electric, the U.S.Weather Bureau, Forsyth Dental School, Boston 
Museum of Fine Arts, Harvard University and Boston University. 

IXirlng the past quarter, there were also 53 Individuals who made brief tours 
of the coaputer installation at different times. Represented by these Individuals were 
the Office of Naval Research, Rsaanual College, Oirtmouth, General Radio Corp., 

Middlebury College, Seoul University, Dean Witter i, Co., W. R. Grace Co., Institute of 
Statistical Mathematics In Japan, Sperry Gyroscope Co., Monsanto Chemical Co., British 
Ministry of Supply. Bell Telephone Mfg. Co., Ramo-Wooldrldge Corp., Rutgers University 
and the U.S. Patent Office. 


SubJ ect 

Description 

Units 

Year 

Instructor 

6.25 

Machine-Aided Analysis 

3-6 

4 

Linvill 

6.535 

Digital Coaputer Coding and 

Logic 

3-6 

G 

Arden 

6.538 

Electronic Coaputational 
Laboratory 

3-1-5 

G 

Versuh 

6.54 

Pulsed-Data Systems 

3-6 

G 

Linvill 

6.567 

Switching Circuits 

3-1-6 

G 

Caldwell 

6.568 

Switching Circuits 

3-1-6 

G 

Caldwell 

2.215 

Methods of Engineering Analysis 

3-9 

0 

Crandall 

15.542 

Management Information Systems 

3-6 

0 

Gregory 

M3 9 

Methods of Applied Mathematics 

3-9 

G 

Hildebrand 

M411 

Numerical Analysis 

3-2-7 

G 

Hildebrand 

M412 

Numerical Analysis 

3-2-7 

u 

Hildebrand 


It Is apparent from the above list that these subjects are predominantly Graduate A 
subjects — which indicates that this subject Is primarily for graduate students. 

However, it Is perfectly possible for undergraduate students In their senior year to 
take these subjects as an elective. 

Machine-Aided Analysis - 6.25 

Subject 6.25 — Machine-Aided Analysis — was offered for the first time during the 
fall term 1953-1954. This subject differs from other coaputational subjects in the 

following respects: 

1. It Is a fourth-year elective subject designed for the undergraduate, 

2. "it Is a survey subject which covers numerical analysis, analog computation, 

and digital computation, 

3. It provides the student with a general knowledge of some of the numerical methods 
necessary for machine computation. 

Admittedly such a survey subject cannot be all Inclusive, and because of the large 
enrollment, it is no longer possible to provide the student with actual physical contact 
with any computer. 

Since 1953 this subject has turned out to be increasingly popular. As a matter of 
fact It Is now offered both fall and spring terms. The attendance is quite high. There 
were'106 rtudents In 6.25 during the Fall Term 1955, and 47 students during the Spring 
Term 1956 This subject offered at the undergraduate level enables students to decide 
whether they are Interested in taking additional subjects in this area, and if they are, 
there are 10 graduate subjects In which they may enroll. 
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Introduction to Digital Computer Coding and Logic - 6.535 

This subject was, as usual, again offered during the Spring Term 1956, and there 
were 97 students enrolled. As a result, It was necessary to divide the subject Into 
two sections taught by D. Arden and F. Hoiwig, respectively. 


Because of the large enrollment In this subject, no attempt was made to make actual 
use of any computer. As a matter of fact, the usual detailed coding for the Whirlwind 
Computer was not presented, and students were not allowed to use the cooputer because 
of the large size of the group. Since the departure of Professor Adams, the content 
of the course has been changed somewhat, and during the spring term considerable 
emphasis has been devoted to a study of the logical design of the Internal system 
aspects of several large computers — Whirlwind I, NORC, XD-1, etc. 

Klectronlc Computational Laboratory - 6.538 

Tlie subject matter offered In 6.538 since 1947 provided the student with two basic 
types of material. Specifically, same 80% of the laboratory work was devoted to a study 
of electronic circuit components by actual laboratory work on the following experiments: 

I. Time Measurement Equipment and Study of Synchroscopes and Oscilloscopes . 

II. Pulse Generating Equipment 

An experiment on blocking oscillators, delay-line pulse generators, 
sharpening and squaring circuits, and ringing circuits. 

III. Bistable Circuits and Flip-Flops 

An experiment covering various flip-flop circuits Including Tbyratron, 
single-tube pentode circuits, Eccle-Jordan circuits, high-speed pentode circuits, 
and transistorized flip-flops. 

IV. Klectronlc Counters and Counting Circuits 

An experiment covering Thyratron counters, binary decade counters, screen- 
coupled counters, biquinary counters, respectively. 

V. Coincidence Circuits and Gate Circuits 

An experiment considering various types of coincidence circuits: diodes, 
trlodes, pentodes, and other multiple coincidence devices. 

VI. Design of Bistable Circuits 

A project-type of experiment In which the student designs (on paper) a 
flip-flop to meet certain operating conditions, and tests his design by actual 
construction of the flip-flop, and testing its performance In the laboratory. 

VII. Magnetic Recording of Pulsed Information 

An experiment on magnetic recording In which the magnetic material head 
recording frequency, pulse amplitude, tape velocity, are varied and the 
resulting performance evaluated. 

VIII. Time Pulse Distributor 

An experiment studying vacuum tubes and diodes as control elements by 
pulse distribution of voltage signals. 

Bquare-HyateroslB-Loop Magnetic Storage Elements 

An experiment on the use of magnetic cores as a storage media, number of 
turns on various windings, current amplitude of recording and sensing turns 
are used as parameters. * ' 


IX. 


Detailed Information regarding the above nine experiments was given to the students 
by mean, of a written set of note. (230 page.) which were written a. reference notes and 


experiment notes, respectively. In addition to the above nine experiments, two laboratory 
experiments were devoted to the use of punched card equipment. Specifically, one on the 
card punch, sorter, and tabulator, and the second on the Card Programmed Calculator, 
respectively. The class lectures wore primarily devoted to a discussion of basic construct¬ 
ion of digital computers, and their use In the solution of scientific computations 
encountered In the solution of simultaneous algebraic equations, inversion of matrices, 
solution of ordinary differential equations, etc. 

Because of the availability of the Type 650 Magnetic Drum Calculator, it was decided 
that the subject content of 6.538 should be revised considerably. Accordingly, during 
the Spring Term 1956 all of the above-mentioned electronic experiments were deleted 
from the subject content and the subject was devoted entirely to the application and ubo 
of digital computing equipment. The students were tuught thw basic principles of digital 
computers by performing experiments on card punches, sorters, accounting machines, tho 
Card Programmed Calculator, and the Type 650 Magnetic Drum Calculator. The principle 
emphasis was actually devoted to a study of the use of the 650 machine in the solution 
of scientific and business data processing problems . 

Students were taught to program the 650 by becoming familiar with the basic 650 
language. After the basic language was understood, they were taught tho following 
interpretive systems: 


1. 

MITSS 

MIT Selective System 


2. 

MITILAC 

FLOATING DECIMAL Interpretive-mnemonic coding 

(1600) 

3. 

NOP I 

Fixed and Floating Interpretive-numeric coding 

(600) 

4. 

SOAP 

Symbolic Optimal Assembly Program 


5. 

FLIMSY 

Floating Interpretive Matrix System (815) 


6. 

BTL 

Fixed and Floating interpretive-numeric coding 

(IOOO) 


It is obvious that the students were not able to become completely familiar with all of 
the above systems. All students Instructed now understand tho 650 language, MITSS, and 
MITILAC. However, only certain students became familiar with each of the systems 
depending upon their particular Interest. 

The following list of homework and experiments were performed during this subject: 

I Problem 1 - Laboratory use of the Type 026 Card Punch and Type 082 Sorter 

The student spends two hours in the laboratory using the card punch and 
sorter to prepare and sort several hundred IBM punched cards, thus insuring 
complete familiarity with these two machines. 

II Problem 2 - Design of Punched Card 8ystem to Handle Business Data Processing 
Work at the Eastern Joint Computer Conference 

This problem teaches the importance of overall flow charts, the design of 
suitable IBM card layout forms, the preparation of accurate instructions for 
machine operating personnel, and the estimation of the personnel and machine 
time required to hrndlo the registration for 2,000 attendees at the conference. 

III. Problem 3 - The Solution of a First Order Differential « aH_»lL on on the Type BSO 
Ma gnetic Drum Calcu lator 
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IV. 

V. 

VI. 

VII. 

VIII. 

IX. 

x. 

XI. 

XII. 

XIII. 
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The student Is given a block diagraw of a 650 program used to solve an 
ordinary differential equation using fourth-order Runge-Kutta formulae. As 
part of the assignment, the student prepares the actual coding of the 650 
program using basic language, prepares a test solution for this program, computes 
a few check points, debugs his program on the 650, and finally carries out the 
solution of this problem on the machine. After completing this particular 
experiment, the student has a very Intimate knowledge of the 650 machine. 


Problem 4 - Solution of an Ordinary Differential Equation on the 6 50 

■me student solves the Identical problem stated lVthF'pT^edlng section 
using the Card Programmed Calculator. This teaches them the difference between 
the internally-stored program features of the 650 and the card-prog rammed 
characteristics of the Card Programmed Calculator. 


ProLlem 5 - Evaluatio n of the Roots of a 98 degree Polynomial 

of 65 ? Pr °?T" U,in * the NOTI 8y8te " prepared rtich evaluates the roots 
of the polynomial using Hitchcock's method. 


Problem 6 - Computation of the Act arlal Life Insurance Tables 

fiin * a hlgh deBr0e °* accuracy ln this computation was needed, the basic 

problem^**** led to perf °™ arithmetic in the solution rflhl. 


rrooiem / 


on 01 “ prices using Modified Elimination Procedure 
A basic 650 program was prepared which automaticallv includeTli'^Dlvotal- 

zerosion® 0 ? t ° Ch " lt ' U8 - correctly handles' matrlce^lcS ha™ 

zero elements on the main diagonal. 

Probtem^8_:^Calculation of Nuclear Transformation, and Radioactive Decay 
in the nth ^ t Pr08 T * BB prepared which computes the number of atoms'present 
in radioactive declv oi 7 8lUUltaneoua transformation by nuclear reaction 
time y 8 V * n t4rget “dide in a given neutron flux at any 


XIV. 


Problem 14 - Iterative Me thod for Solving Simul taneous equations on the 650 
A 650 Program was prepared which evaluates the convergence rate, 

accuracy and solution time using Craig's method, Conjugate Gradient method, 
and Lanczos's method, respectively. 


XV. Problem 15 - 650 Solution of the Vertica l Resistivity Problem 

A MITSS program was prepared for the solution of the Integral equation 
arising In geophysical prospecting. 

XVI. Problem 16 - Preliminary Analysis of an Inventory Control System 

A non- machine study of the many problems arising in a systems study of 
the overall Inventory control problem. Detailed description Is given of the 
nature cf transactions involved, considerations Involving volume of data, rate 
of Information flow, need for rearrangement of computational procedure, study 
of the offects of errors, breakdown, and Interruptions, and finally a consid¬ 
eration of problem areas of interest to management to provide better overall 
control. 

XVII. Problem 17 - Solution of a Payroll Problem on the 650 Calculator 

A basic 650 program was prepared for a hypothetical payroll Involving 
200 employees. The study Included a consideration of piece-work as well ar 
the regular payroll in which employees switch back and forth several times 
each day. The usual Items such as, regular pay, overtime pay, bond deduction, 
hospitalization, Income tax, union dues, credit union, Community Chest, and 
other miscellaneous deductions, were included. 


XVIII, Problem 18 - Hitchcock's Method for Bvaluatlng the Roots of a Polynomial 

A basic 650 program was prepared for evaluating the _ roots~of~a~poiynomlal. 
Particular emphasis was given to a consideration of polynomials which cannot 
by handled by this method. 


TihiuTT interes t Calculations for Savings Banks 

• lx difSfSt "ttoST WBS Pr ° ParSd WhlCh “'optically handTes - any of the 
^different method, for performing Interest calculations arising in the savings 

Pr ° b 7blslc' B1Varlate Fr " nUfln ^ ^tributlons on the 650 

distribution by corre,pon"ng P 77^7ls7“7.7 blvaria t« frequency 

also computes certain statistic . (tor various parameters), and 

An effective method of storage AllMaMon 7 above - Bontion «<l distribution, 

for obtaining the desired argument, and t * OBplo >’ od employs a neat trick 

of values. 1 d functlon » in an Internally-stored table 

Pr° bl e»bll - Study of Optimum Programming on the 650 Calculator 

semi-optimum data^r!!^,”" tbe^T! Wt ’7 * tUdl8d variou « "ethods for effecting 
selection using T^S^lS, 7"" 

Problem 12 - Solvent Extraction of lletal. 

of zircon~~and hLn£T“ "* ^ouT'a^or^c^a^.^distribution coefficients 
Pr ° bl :W 7 UtiOD of the "■»»* Plow Equation 

equatlon~de^criblng 0 the 111 heat floITproblea usl 80lUtl0n 0t a partlal differential 
conditions. Problem using several different boundary 


XIX. Problem 19 - Study of Series-Paralle l Switchi ng F unctions 

A basic 650 program was prepared for evaluating the switching functions 
of six or fewer variables. The program employed is a modification of the 
UcCluskey ^iine reduction procedure by means of which the preliminary lapll- 
cants and minimum sum forms of the switching functions were obtained. The 
solutions of 86 different problems were performed during the course of this 
study. 

XX. Problem 20 - Solution of Spheroidal Antennae Equation Using the 650 

A NOPI program was prepared for the solution of this equation using 
Milne's predictor-corrector formulae. 

XXI. Problem 21 - 650 Solution of the Production Scheduli ng Probl em 

A MITSS program was prepared for the solution of this problem which 
determined production curves which minimize costs for a single arbitrary 
sales forecast. Graphical results were obtained which related the sensitivity 
of the solutions for various system parameters using a quadratic production 
cost function. 

XXII. Problem 22 - Solution of the Inventory Control Problem on the 650 Calculator 
A basic 650 program was prepared for a simplified inventory control 
problem. Specifically, the program was designed to handle any number of stock- 
numbered items up to 34-per-day items. 'Hie program was limited to control and 
reporting of customer products only. No attempt was made to control the com¬ 
ponent parts for the items. Similarly, no attempt was made to report the raw 
material requirements. 
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XXIII. Problem 23 - Study of Brokerage Accounting 

A basic 650 program wh prepared for calculating the items appearing In 
the transaction statement of a customer's brokerage account. The following 
items were Included in this study: 1) Commercial postage and Insurance, 

2) Federal transfer taxes, 3) State transfer taxes, 4) S.B.C. fees. 

Naturally a program of this type cannot be all inclusive and a number 
of simplifying assumptions had to be made to effect the solution of the 
program in a reasonable amount of time. 

A consideration of the above list reveals that a wide variety of problems were 
considered. It must be emphasized that the latter group of subjects were performed on 
a term paper" basis. That is, these studies were performed on an individual basis, and 
elaborate reports were prepared describing the problem statement, the actual 650 program, 
the results obtained, and a summary and evaluation of the paper. Although the term 
paper was expected to require 30-40 hours of home and laboratory work, the actual time 
put in by the students varied enormously. For example, this time ranged from 30 hours 
to 360 hours, respectively. In the latter case, students intend to use this topic as 
a basis for their master's theses during succeeding terms. 

It is apparent from the above course description that an attempt was made to cover 
both scientific computing and electronic data procoasing in the subject. This can 
indeed be done, particularly when individual term papers are assigned. However, as time 
goes on, it may be advisable to separate this subject into two distinct subjects, one 
dealing solely with scientific computing and the other with business data processing. 

This subject was taught by Dr. F. M. Verzuh, Director, Office of Statistical Services. 
There were 26 students enrolled in Subject 6.536 during the Spring Term 1956 of which 
there were 12 Electrical Engineering students, 4 Chemical Engineering students, one 
Geophysics student, 3 Economics students, 2 Business Administration students, and 4 
Mathematics students, respectively. This was the first time students were able to use 
the 650 Calculator as part of their laboratory work. Judging by the large number of 
this group — 12 students — who are now using the machine for the solution of their 
thosis problems, it appears that the subject was highly successful. 

In addition to the solution of scientific computing and business data processing 
problems, a number of classroom lectures were devoted to the consideration and evaluation 
of the relative merits of the many commercially-available electronic data processing 
systems. Specifically, this was achieved by a study of 16 representative intermediate- 
size drum computers, and seven of the large-scale systems — B17.MAC, Datamatic 1000, 

ERA 1103, IBM 704, IBM 705, Raycoo, and UNIVAC, respectively. Since many of the students 
were entering various large business, industrial, and government concerns, it was 
important that they have some knowledge of some of the relative merits of the individual 
systems. Finally, several lectures were devoted to the IBM 704 Calculator which will be 
Installed In the M.I.T. Computation Center in the near future. 


Subject 

Nome 

No. 

Fall 

Students 

spring 

Total 

2.215 

Methods of Engineering Anulyflis 

18 

0 

18 

6.25 

Machine-Aided Analysis 

106 

47 

153 

6.535 

Digital Computer Coding and Logic 

0 

97 

97 

6.538 

Electronic Computational Lab. 

0 

26 

26 

6.54 

Pulsed-Data Systems 

0 

60 

60 

6.567 

Switching Circuits 

48 

0 

48 

6.568 

Switching Circuits 

0 

30 

30 

15,542 

Management Information Systems 

0 

22 

22 

M39 

Methods of Applied Mathematics 

130 

123 

253 

M411 

Numerical Analysis 

45 

0 

45 

M4I2 

Numerical Analysis 

0 

16 

16 


Total 

347 

421 

768 


It is apparent from the above Hat that the total attendance in those 11 subjects was 
768 students. This does not mean, however, that there were 768 individuals involved 
Since many of the students took a number of the above subjects, thus the actual number 


V. M. Verauh 

Office of Statistical Services 


Su mmary o f St udent Enrollment Statistics 

A detailed description of the course content of the subjects listed in the intro¬ 
duction is given in the M.I.T. Catalogue. 

An indication of the interest in machine computation may be obtained by considering 
the following enrollment in each of these subjects during the academic year 1955-1956: 


70 


71 






APPROVED FOR 


PERSONNEL OF THE PROJECTS 


MACHINB METHODS OF COMPUTATION AND NUMERICAL ANALYSIS 

Faculty Supervisors : 

Philip M. Morse, Chairman 
Samuel H. Caldwell 
Herman Feshbach 
Jay W. Forrester 
James B. Reswlck 
Chla Chlao Lin 

Research Associate : 

Bayard Rankin 

Research Assistants : 

Joseph Hershenov 
M. Douglas Mcllroy 
James W, Schleslnger 
Hurbertus J. Weinitschke 
Fernando J. Corbato 
Harvey Fields 
Zcltan Fried 
Raymond F. Stora 
Arnold Tubis 
Aaron Temkin 
Jack L. Uretsky 
Jukka A. Lehtinen 
Joseph B. Walsh 
Marius Troost 
Norman F. Ness 
Ralph 0. Gray 


PROJECT WHIRLWIND 


Staff Members of the Scientific 
Digital Computer Laboratory: 


and Engineering 


Frank M. Versuh, Head 
Dean N. Arden 
Sheldon F. Best (Abs.) 
Frank C. Helwlg 
Leonard Roberts 
Jack Roseman 
Arnold Siegel 
Mirray Watkins 
Munroe R. Weinstein 
Alexander Zabludowsky 


Physics 

Electrical Engineering 
Physics 

Electrical Engineering 
Mechanical Engineering 
Mathematics 


Mathematics 


Mathematics 

Mathematics 

Mathematics 

Mathematics 

Physics 

Physics 

Physics 

Physics 

Physics 

Physics 

Physics 

Mechanical Engineering 
Mechanical Engineering 
Chemical Engineering 
Geology 

Civil Engineering 


Computations Group at the 
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